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Abstract The energy and particle fluxes emitted by an accelerated two level atom 
are analysed in detail. It is shown both perturbatively and non perturbatively that the 
^ ■ total number of emitted photons is equal to the number of transitions characterizing 

thermal equilibrium thereby confirming that each internal transition is accompanied by 
the emission of a Minkowski quantum. The mean fluxes are then decomposed according 
to the final state of the atom and the notion of conditional flux is introduced. This notion 
is generalized so as to study the energy content of the vacuum fluctuations that induce 
the transitions of the accelerated atom. The physical relevance of these conditional 
fluxes is displayed and contact is made with the formalism of Aharonov et al. The 
same decomposition is then applied to isolate, in the context of black hole radiation, 
the energy content of the particular vacuum fluctuations which are converted into on 
mass shell quanta. It is shown that initially these fluctuations are located around the 
light like geodesic that shall generate the horizon and have exponentially large energy 
densities. Upon exiting from the star they break up into two pieces. The external one is 
red shifted and becomes an on mass shell quantum, the other, its "partner", ends up in 
the singularity. We avail ourselves of this analysis to study back reaction effects to the 
production of a single quantum. 
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1 Introduction 



The history of Hawking radiation and of uniformly accelerated detectors in Minkowski 
vacuum have to some extent evolved in parallel. Indeed Hawking's discovery |TJ] that 
black holes emit particles at temperature 1/8ttM (M being the mass of the black hole) 
closely preceded and inspired Unruh's discovery^ that a uniformly accelerated detector 
will thermalize at a temperature a/2ir (a being the acceleration of the detector). The 
global structure of the geometry (the presence of horizons), the singular behavior of the 
modes and the mechanism of particle emission (excitation of the detector) are in close 
analogy in both problems. 

In this article we exploit this isomorphism to obtain a quantum description of the 
fluctuations of the energy density in both problems. To this end, we first describe the 
energy content of the field configurations correlated to excitations of an accelerated 
system. We then apply the same analysis to obtain the energy distribution of field 
configurations (the vacuum fluctuations) which get converted into Hawking quanta. The 
main point of this analysis is that it goes beyond the description of the mean (averaged) 
energy of the field (see ref. || for a review): we compute the energy of the fluctuations 
around the mean and evaluate specific gravitational back reaction effects induced by 
these vacuum fluctuations. 

In Part 2, following the work of references P|H]-[ffO|, we analyze the mean energy 
density emitted by a uniformly accelerated detector. We emphasize the difference be- 
tween global quantities (the total energy, the number of quanta) and local quantities (the 
energy density) and we exhibit the essential role of transients in ensuring that global 
properties are respected (such as the positivity of the total energy). Our results gener- 
alize previous perturbative||10|| and non perturbative|7]] analyses. In particular, we show 
that the total mean number of Minkowski photons emitted is equal to the mean number 
of internal transitions of the two level atom which have occurred during the interacting 
period. 

In Part 3, we decompose the mean flux analyzed in Part 2. It is written as a sum of 
two terms 0] || [Kj. The first term is the energy density if the atom has made a transition 
and the second term is the energy density has not made a transition. From these energy 
densities we introduce the central notion of conditional energy emitted. This notion 
is then generalized so as to investigate others correlations. We decompose the mean 
vacuum energy (which is zero) into the energy densities correlated to a future transition 
(or to the absence of a future transition) of the two level atom. These conditional 
energy densities describe the vacuum fluctuations which shall give rise to transitions of 
the atom. 

This description follows the treatment of ref. |TI| wherein the mean current density 



carried by created pairs in an external electric field was decomposed into conditional 
current densities so as to describe the current carried by a specific pair of quanta. 

As this approach to isolate certain vacuum fluctuations based on conditional values of 
operators is rather new and leads to peculiar results such as complex energy densities we 
shall dwell on it somewhat. In Section [0]we show how the non diagonal matrix elements 
which describe these vacuum fluctuations can be obtained by decomposing the mean 
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value of the energy density. Their physical relevance are then revealed by introducing 
an additional quantum system coupled to T^ v . Indeed, these complex densities have a 
natural interpretation in Quantum Mechanics only. In perturbation theory, one finds 
that the modification of the wave function of the additional system is controlled by these 
matrix elements. (In particular, in the black hole situation, the first order response 
of the metric to the creation of a particular Hawking photon is controlled by these 
matrix elements.) This was discussed by Aharonov and collaborators |T2[ who considered 
the additional system to be a measuring device as in the von Newmann approach to 
measurement theory fT3H. For the sake of completeness, at the end of this article, in Part 
5, we show how the approach of Aharonov et al. can be adapted and generalized to the 
present cases where we are more concerned by backreaction effects than measurability 
problems. 

The analysis of vacuum fluctuations is applied to black hole radiation in Part 4. We 
obtain the energy density of the vacuum fluctuations which are converted by the time 
dependent geometry of a collapsing star into Hawking photons. The properties of these 
energy densities give rise to specific quantum gravitational back reaction effects which are 
not present in the mean theory. At the present time, the only mathematically consistent 
treatment of the backreaction is the semiclassical theory (see refs. [|T1J[|T^(TB|]). In this 
approximation, the external field remains purely classical and only the mean value of the 
matter current operator acts on it as a source. The quantum properties of the matter, 
i.e. its fluctuations and correlations, are completely ignored. But it is the fluctuations 
which are problematic. Indeed it has been stressed by t'Hooft[17| and Jacobson [18| that 
Hawking's derivation of black hole radiation is no longer valid as soon as gravitational 
interactions are taken into account because it makes appeal to the structure of the 
vacuum on exponentially small scales. In the present work this is seen in particularly 
vivid fashion. Indeed we shall show that if a particle is emitted by the black hole a 
time u after collapse with asymptotic frequency A, then there was a vacuum fluctuation 
inside the star of energy density 0(u 2 ) (where u = Ae U//4M ) located on a distance scale 
of order uj^ 1 . After a time u = 0(4M In M) for a typical A = 0(M _1 ), u is greater than 
a Planck frequency and the free field theory of Hawking is unjustified. This aspect of 
the fluctuations has been presented together with F. Englert in ref. [|i~9|j . 

There have been two attitudes in the literature to confront this situation. The first 
is to try to guess what is the physics at the Planck scale near the horizon and how the 
Hawking radiation emerges therefrom |[21||-||25||. The second has been to use Einstein 
equations to investigate how back reaction effects modify the production of Hawking 
photons p6|[|27 |. The present article places itself in this latter vein. We show at the end 
of Section [O] how our results can be applied to study some simple back reaction effects. 
We first evaluate the change in probability of finding a specific Hawking photon due to 
a modification of the background metric. We then show how the creation of a particular 
photon modifies the probability of finding subsequent photons. In both cases, one sees 
explicitly that the change in probability is entirely controlled by the conditional value 
of the energy momentum tensor. The more difficult problem of the self interaction of a 
Hawking photons with itself as it is created necessitates the analysis of loop corrections 
and we hope to report on it in a subsequent paper. 



3 



2 The Energy Emitted by an Accelerated Atom 



2.1 Introduction 

It is now well known that a uniformly accelerated two level atom thermalizes in Minkowski 
vacuum at temperature a/27r@. But it is much more complicate to obtain a complete 
description of the fluxes emitted by this thermalized atom. 

As first pointed out by Grove||, when thermal equilibrium is reached there is no 
net emission of energy. His argument is the following. The accelerator feels the effect 
of a thermal bath. So first consider the inertial two level atom in thermal equilibrium. 
The time independence of the Hamiltonian and the stationarity of the state of the atom 
(in the thermodynamic sense) guarantee no mean flux, since each absorbed photon 
is re-emitted with the same energy. This argument is immediately applicable to the 
accelerator since the fact that a = constant implies that his physics is translationally 
invariant in his proper time r (i.e. invariance under boosts). Since Minkowski vacuum 
is also an eigenstate of the boost operator, the total Rindler energy is conserved. 

A detailed picture of the steady state emerges from the following consideration. Focus 
on the ground state of the accelerator which excites by absorbing a Rindler quanta (a 
rindleron) coming in from its left. Then the field configurations to its right is depleted of 
this rindleron. Since this rindleron carried positive energy, its removal can be described 
as the emission of negative energy to the right. In equilibrium there is also to be 
considered the process of disexcitation corresponding to the emission of positive energy 
to the right. The Einstein relation guarantees that the two cancel. This implies no net 
energy flux. 

However the accelerated atom is in Minkowski vacuum, hence all perturbations of 
the radiation state lead to the production of Minkowski quanta. Then, how do we 
reconcile the (certainly) positive energy of these produced quanta with the absence of 
radiated energy in thermal equilibrium? The answer lies in a global treatment of the 
radiation field which also takes into account the transients due to switching on and off 
the detector (or equivalently the transients which occur when the detector passes from 
an inertial trajectory to the accelerated one). 

That transients may have a global content which depends on the whole history also 
occurs in the problem of the classical electromagnetic field emitted by a uniformly ac- 
celerated charge in 3 dimensional Minkowski space (see ref. P%|]). We point out that 
this situation is particular to uniformly accelerated systems and makes explicit appeal to 
the exponential Doppler shift between accelerated and inertial reference frames. Indeed 
(anticipating some notations introduced in the next sections) one can express the total 
Minkowski energy emitted (in V-modes) as 

r+oa r+oo (In, 

E M = / dV(T vv ) = / dv(T vv )— = Qe" 074 - GjeT^ (1) 

JO J— co dv 

where the first integral is limited to the domain V = t + z > because the accelerator 
remains in the quadrant z > \t\. In the second equation we have used the Rindler 
coordinate v related to V by dv/dV = e~ av ; and in the last equation we have made 
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appeal to Grove's theorem which states that the integrand vanishes everywhere except 
at the endpoints (= r») and Vf (= r/); Cj and C/ depend on the exact form of the 
transients at Tij. 

Equation (HD can also be written as an integral over the rates of absorption and of 
emission of Rindler photons taking into account that a transition at time r is accom- 
panied by the emission of a Doppler shifted Minkowski photon of frequency ui{r) = 
mdv/dV = me~ aT where m is the resonant frequency of the atom. Hence the energy 
emitted is 

Em = dr Ru{t) 

•>Ti 

= ^(e- ar > - e- aT f) + (C' ie - an - C' f e- ar f) (2) 

Here R is the number of transitions per unit proper time. The constants C[ and Cj 
depend on the detailed way the interaction is turned on and off at nj. So the integral 
eq. (§) is of the same form as eq. (jH). These two expressions are compatible because 
all the photons in eq. (|2|) interfere in such a way that their energy is only found at the 
edges of the interaction period as in eq. ([!]). 

This part is organized as follow. Section 2.2 is devoted to recalling the main properties 
of Rindler quantization. In Section 2.3 we present the model of the accelerated two level 
atom. In Section 2.4 we obtain formal expressions for the energy emitted. We decompose 
the flux according to the final state of the atom thereby introducing the central notion 
of conditional energy emitted. The properties of these conditional fluxes as well as their 
physical meaning shall be displayed in Part 3. As a warm up, we present in Sections 
2.5 - 2.7 the various properties of the mean flux insisting on the role of the transients 
in guaranteeing that global properties are respected. The mean fluxes as the atom 
thermalizes are discussed in Section 2.5 and then in thermal equilibrium in Section 2.6. 
Finally in Section 2.7 we use the exact solvable model of ref. to extend the previous 
perturbative results to all order in g. 



2.2 The Rindler Quantization in 1+1 Dimensions 

In this section, we review the relevant properties of the Rindler quantization of the 
scalar field in Minkowski space time in 1+1 dimensions. The conformal invariance of the 
massless scalar field is best exploited by using the light like coordinates U, V defined by 

U = t-z 

V = t + z (3) 

Whereupon the Klein-Gordon equation takes the form dxjdycj) = and any solution can 
be written as 

mv) = 4>(U) + </>(V) (4) 

From now on we shall drop the right moving piece and consider the "V" term only. It 
is obvious that all conclusions shall be equally valid for the right movers. 
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The second quantized field can be decomposed into the orthonormal complete basis 
of Minkowski modes 

0(V) = J™ du (auVu{V) + (5) 

-iOJV 

MV) = -7= (6) 

Minkowski vacuum |0m > is the state annihilated by all the a^'s. The propagator 
(Wightman function) in Minkowski vacuum is 

G+(y,V) = <o M \(f>(v)(f>(v , )\o M > 

roc 1 

= ^ d^iy^V) = -—\og{V -V -ie) (7) 
The (normal ordered) hamiltonian of the field is (for left movers) 

r-\-oc poo 

Hm = \ dVTw = I dujuj(aljauj) (8) 



where 

Tyy = <9y0<9y0 (9) 

Therefore Minkowski vacuum is the ground state of the hamiltonian Hm'- Hm\^m > = 0. 

The uniformly accelerated observer will be taken to be in the right (R) Rindler 
quadrant U < 0, V > 0. In this quadrant one defines Rindler coordinates p, r by 

j t = p sinhar , . 

1 x — p coshar 

The accelerated observer follows the trajectory p—l/a (where a is its acceleration) and 
its proper time is r (see figure 1). In this quadrant one introduces also the light like 
Rindler coordinates u, v defined by 

u = t — a^ 1 \nap 

v = r + a~ 1 \iaap (11) 
These are related to the light like Minkowski coordinates U, V by 

f U = -a- 1 e' au 



V = a 



(12) 



The coordinates eq. flTU| ) may be extended to the left (L) Rindler quadrant by the 
analytic continuation r — >• r ± in /a. One introduces in L the light like coordinates 
u L ,v L given by 

J U = a~ l e aUL , , 

V = -a- 1 e- a ^ [L6) 
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The natural basis of quantization the uniformly accelerated observer would choose 
is the Rindler basis which consists of plane waves in the variables u, v (Rindler modes). 
This is because u, v are related to its proper time r as U, V are related to the Minkowski 
time t. The Rindler v modes are thus given, in strict analogy with eq. (|j), by 

g— iXv 

= 75a (14) 

Since the <px,R constitute a complete set in R (V > 0) only, they cannot be related 
to the Minkowski basis by a unitary transformation. One must also introduce Rindler 
modes living in the left quadrant. But one finds that the Bogoljubov transformation 
relating the Minkowski modes to the Rindler modes is singular at V = |2£| and care 
must be taken to define it as a limit if the Minkowski properties of the theory are to 
be satisfied. To this end it is useful to first introduce an alternative basis of positive 
frequency Minkowski modes, eigenmodes of iaVdy (= id v for V > 0), and defined for 
all V 101: 



[a(e + iV)]~ iX/a 



= e^o 



^( e 7rA/a _ e - 7rA /a)47TA 

e nX/2a Q(V)(aV)-" lX l a ^ e -^A/2a q(_ (_ a y}-i\/a 



where 



irX/a _ e -nX/a\ ,/47r|A| /I gTrA/a _ g -7rA/a I J47r|A 



i\/a 



(15) 



an \ 1 (u 



^ T(i\/ a) \l XsmhirX/ a) ^2nacu la ) 



The first factor in j\,oj is a pure phase introduced for convenience. The factor e~ ue 
is the crux of the construction. It defines the integral eq. (0), regularizes the modes 
tpx.Niiy) at V = and ensures the correct Minkowski properties of the theory. For 
instance it gives the correct pole prescription at V — V of the propagator eq. (^) when 
expressed in terms of the modes <fx,M as G + (V, V) = J_™ d\(px,M(V)(p* x M (V'). The 
limit e —> is to be taken at the end of all calculations. 

The annihilation and creation operators corresponding to the modes <f\ t M are a x ,M 
and at 



b X,M- 

f +00 



<P(V) = / dX Lx,mVx,m{V) + 4,aM,m(V)) (17) 

The right and left Rindler modes <P\,r(V) and (Px,l{V) can now be defined by the 
linear unitary transformation 

UxM = ^n + ^l,L A>Q (lg) 
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with 



ir\/2a —ir\/2a 

a A = =_ , Px= , x/ =fr and a\ - $ = 1 (19) 

^/g7rA/a g— TTA/a «/ g7rA/a g — nA/a 



In the limit e — ► the Rindler modes take the familiar form (see eq. (0)) 



p AiL (V) = ^(-V)(-a1/) iA / a /V4^rA = e-^/v^rA 



(20) 



For finite e they differ from these limiting forms only when V < e. The Rindler de- 
struction and creation operators ax,R,a\ R and a\ t L,a\ L associated to these modes are 
related to the Minkowski operators a x ,M by the following Bogoljubov transformation 

{ a \,R = «A«A,Af + P\ a -X,M ^21) 
o-\,l = a A a_A,Af + P\ a X,M 

(by virtue of the orthonormal character of the two sets of modes (px,R, <p x ,L and <p x ,M as 
well as a\ — (3\ = 1). One immediately deduces that the mean number of Rindler quanta 
present in Minkowski vacuum is given by the Bose Einstein distribution 

< M \4 jR a x ,, R \0 M > = — L- ^(A - A') (22) 

It is useful to introduce the generator of boosts H R since it generates translations in r 
and is therefore the Hamiltonian for the accelerated observer. 

/+oo 
dVaVT vv 
-oo 

r+oo r+oo 

dvT vv - / dvrX 

-oo J —oo 

r+oo 



LJ-v L v L 







d\X(a\j i ax j R — a\ L a\,L) 
d\X(a^ M ax t M ~ a -A,M a -A,A/) (23) 



One sees that quanta in L carry negative "Rindler energy". This is because r (defined 
by arctanh(t/z)) goes backwards in time in L (dr/dt < in L). Furthermore since 
Minkowski vacuum is invariant under boosts it is annihilated by the Rindler energy 
operator: 

H R \0 M > = (24) 
Using the above Bogoljubov transformation it is easy to show that 

\0m > = II— e~k<^\,n\Q RL > (25) 



A 



a>x 



where \0rl > = |0_r>®|0l> is Rindler vacuum in both the right (R) and left(L) 
quadrants. Thus, from eq. (p3|) , one sees that the pairs of Rindler quanta (whose mean 
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number is given by eq. (|22|)) present in Minkowski vacuum carry zero Rindler energy 
One sees also that upon tracing over the left quanta the reduced density matrix in the 
right Rindler quadrant is an exact thermal distribution of right rindlerons [[|. 

Eq. ( p5|) will be useful in the next chapters since it provides an easy way of calculating 
the probability (| (i[>r\0m) | 2 ) of finding in Minkowski vacuum a state ( \ipn >) containing 
a given number of Rindler quanta. 



2.3 The Uniformly Accelerated Two Level Atom 



The situation we consider is a uniformly accelerated two level atom coupled to the 
massless field cf) introduced in the previous section. The trajectory of the two level atom 
is given by eq. (|I0|) with p = a" 1 : 



t a (r) = a 1 sinhar 
VJt) = 



a- l e aT 



UJt) 



a 1 coshar 



(26) 



The time integral of the interaction hamiltonian is 



/ 



dtdx H int (t, x) 



gm 



dr 



(f{r)e~ mT A + r(r)e im Mt) 0(t a (r), x a {r))] (27) 



where g is a dimensionless coupling constant that shall be taken for simplicity small 
enough that second order perturbation theory be valid, m is the difference of energy 
between the ground and the excited state of the atom, A is the lowering operator that 
induces a transition from the excited state to the ground state of the atom and f(r) is a 
dimensionless function that governs when and how the interaction is turned on and off. 
The factor m on the r.h.s. of eq. ( p7|) is introduced for dimensional reasons. In addition, 
we shall assume that the V-pa.it of the field only is coupled to the atom. This is a 
legitimate truncation owing to the chiral character of the field in 1 + 1 dimensions. For 
simplicity of notation it is convenient to rewrite eq. ( ^7|) as 



dtdx H int (t, x) 



gm 



where 



dre +lmT f*(r)(t)(T) 
We shall be most interested in the situation where /(r 



(28) 



(29) 



Ti 



< r < Tf and /(t) tends to zero outside this interval. 



1 inside a long interval 
Then in the limit Tf — Ti = 



T — > oo with g T finite (i.e. in the Golden Rule limit) the concept of a transition rate 
emerges and is due to the resonance of the Minkowski vacuum fluctuations with the 



fixed Rindler frequency m [30]. The fundamental reason why we need to work with an 
explicit switch function is that we want finite energy momentum densities everywhere 
including the horizon. We shall show that this is possible only for sufficiently rapidly 
decreasing /(r) (which amounts to deal only with wave packets of Rindler modes which 
are well defined in the ultraviolet -see eqs. ([15|, [IB])). 
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When e~ imT /(r) contains no negative frequency in its Fourier transform with respect 
to r, eq. (p8|) defines a Lee model: were it inertial it would only respond to the presence 
of Minkowski particles. For a Lee model /(r) must necessarily decrease less rapidly 
then an exponential when r — > ±00. Alas this condition is too strong and will lead 
to singularities on the horizon in the uniformly accelerating situation. Hence we shall 
be obliged to consider non Lee models which can spontaneously excite. However by 
choosing /(r) such that the negative frequency part of e~ imT f(r) is exponentially small 
the spontaneous excitations are exponentially . Moreover the spontaneous excitations 
occurs only at switch on and switch off transitory periods but do not occur during the 
steady regime when /(r) = 1, i.e. these excitations do not contribute to rates. 

Let us consider first the situation in which both the atom and the field are initially 
in their ground state. The state \ip- > at t = —00 is thus 

\ij)_(t = -00) > = |0 M >|- > (30) 

where |— > (|+ >) designates the ground (excited) state of the atom. At t — +00, when 
the interaction has been switched off, the state can again be expressed in terms of the 
uninteracting states and it is given, to order g 2 , by 

|^_(t = +oo)> = 



where 

P = l 

2 

and where e(r 2 — tl) 

in order to isolate the steady regime part of the interaction. The latter is given by the 
term proportional to ^In^rn whereas the D term is concerned with the transitory periods 
associated with the switch on and off. Indeed to order g 2 the T> term does not contribute 
to the energy density emitted in the steady state regime (i.e. its energy density scales 
like g 2 /T rather than like mg 2 ). Furthermore it carries no Minkowski nor Rindler energy. 
This is proven in the Appendix. We shall therefore drop this term in the sequel. To get a 
flavor of this, the reader can already verify that T> does not contribute to the probability 
to remain in the ground state. Only the term proportional to <$ m <\>m, does so. 

The probability P e for the two level atom to get excited is, in second order pertur- 
bation theory 

P e = g 2 m 2 < M \<f>l<f> m \0 M > (33) 
10 



|0 M >|- > 

/+00 
rfrr(r)e +mr 0(r)|O M >|+ > 
-00 

/ + OO FT 
dr / dT'f{T)e- imT (j){T)f*{T')e +imT <f>{r')\0 M >|- > 
-00 J— 00 

|0m >h > - igm4> m \0 M >|+ > - -g 2 m 2 (pernio M >|- > 
-g 2 m 2 V\0 M >\-> (31) 

J dr 2 J d ri e{r 2 - r^e^ /(r 2 )0(r 2 )e + ^ /*(r 1 )0(r 1 ) (32) 
= 9{j2 — Ti) — 9{ji — r 2 ). We have split the g 2 term in two pieces 



Had we coupled the atom to both the U and V parts of 0, the probability would have 
been twice P e . When /(r) is equal to 1 between n and tj and r/ — Tj = T — > +00 while 
g 2 mT remains finite, the operator 0j„0 m appearing in eq. ( j33|) becomes the counting 
operator for rindlerons of energy m (= a) mR a m ^R) multiplied by n/m. A direct golden 
rule calculation then shows that the probability for the uniformly accelerated atom to 
get excited is 

P e = (l/2)g 2 mTN m (34) 

where N m = i/( e 27rm /« _ 1) i s the mean number of Rindler quanta present in Minkowski 
vacuum, see eq. (0). This proves that the atom maintained on an accelerated trajectory 
reacts to the mean number of Rindler quanta as the same atom, put on an inertial 
trajectory, would have reacted to the mean number of Minkowski quanta. 

If the initial state is the product of Minkowski vacuum and the excited state for the 
atom 

\ip+(t = -00) > = \0 M >|+ > (35) 

then \ip + > at t = +00 is 

\i[) + (t = +00) > = \0 M >\+ > - igmcj^^M >|- > - -^m 2 </> m </4|0 M >|+ > 

+g 2 m 2 V\0 M >\+ > (36) 

where the operator T> which appears is the same as in eq. ([U]). The probability to be 
found in the ground state at t — +00 (i.e. the probability of disexcitation) in the golden 
rule limit is 

P d = g 2 m 2 < O A/ |0 m 0tjO A/ > = (l/2)g 2 mT(N rn + 1) (37) 

Hence, at equilibrium, by Einstein's famous argument, the probabilities P + , P_ to be 
in the excited or ground states are given by 

P P N 

t± = tJL= iVm = e ~ 2 ™/ a (38) 
P_ P d N m + 1 y ' 

that is a thermal distribution at temperature T = a/27r. 
2.4 The Mean Fluxes to order g 2 

We investigate now the properties of the flux emitted by the accelerated atom and its 
relations with the transition probabilities P e and Pd- Since the field is massless, the 
conservation law reads 

duT vv = (39) 

Thus since the interaction, eq. (PTD, occurs only on the accelerated trajectory (eq. 
one has in the past of this line, or to the right of it by virtue of eq. (|39"D, 



(Tw{V,U<U a (T))) = (40) 
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The mean value means that the initial state can be any combination of > and >. 
(When the initial state does contain Minkowski quanta, the mean value eq. (|40|) should 
by understood as the modification of the mean induced by the coupling eq. (]28|).) 

We consider first the situation where the atom is in its ground state at t = — oo 
and the field in Minkowski vacuum (i.e. the state of the system is >)• After the 
interaction is switched off, the mean flux emitted on X + (i.e. on U = +00) or equivalently 
to the left of the trajectory by virtue of eq. fl5§D is, to order g 2 , 



(T VV {V))^_ = <^{t = +oo)\T vv {V,U>U a {r))\^{t = +oo)> 

= g 2 m 2 < Q M \4>m T VV<l>m\^M > - 9 2 m 2 Re < ^M^VV^m^m^M > (41) 



where we have used eq. ( plf) and dropped the contribution of the T> term. The physical 
meaning of the two terms on the r.h.s. of eq. ([II]) was first discussed in M. Before 



discussing it we first rewrite eq. ([41]) as 

(Tvv(V))tp_ = P e (T vv ) e + P g (T vv ) g (42) 

where we have defined 

(T vv ) e = g 2 m 2 <Q M \(t) ] m T VV (t) m \QM>/Pe 

(Tvv) g = -9 2 m 2 Re [< ^ M \Tvv4>l<Pm\^M >] /P g (43) 

where P e and P g are the probabilities to find the atom in the excited or ground state at 
t = +00. P e is given in eq. (]34]) and P g = 1 — P e . 

The interpretation of the two quantities {T vv ) e and (T vv ) g is clear when one recalls 
their origin. (Tw)e comes from the square of the second term of eq. (|3lf) (linear 
in g) whereas (Tw)g comes from an interference between the first term of eq. (|5T| ) 
(unperturbed) and the third term (in which the interaction has acted twice). Hence 
(T vv ) e is the energy emitted if the atom is found excited at t = +00 and (T vv ) g is 
the energy emitted if the atom is found in the ground state at t = +00. These fluxes 
have been normalized so as to express the r.h.s. of eq. ([42]) as the probability of finding 
the atom in a final state times the energy emitted if that final state is realized. Thus 
(Tw)g and {Tyv)e are the conditional "mean" energy emitted. (The word "mean" is 
understood here in its quantum sense, i.e. as the average over repeated realizations of 
the same situation: the same initial state \ip_ > and the same final state of the atom 
-see Part 5 for further comments on this point). 

Similarly, when the initial state of the system is >, the mean energy emitted is 

(Tvv(V)U + = <Mt = +oo)\T vv (V,U>U a (r)M + (t = +oo)> 

= g 2 m 2 < $ M \<p m T V v4>U®M > ~ g 2 ™ 2 Re [< O M |lW0 TO 4jO M >1(44) 



where we have used eq. fl36|) and dropped the V term as well. As in eq. (fl2"D, we rewrite 
this flux as 

(Tw)^ = P d (T vv ) d + P h (T vv ) h (45) 
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where Pd is the disexcitation probability given in eq. fl3?D and where is the probability 
to be found in the excited state at t = +00, hence P^ = 1 — Pd- The conditional fluxes 
(T vv ) d and (T vv ) h are given by 

(Tw)d = g 2 m 2 < Q M \(frrnTvV<Pln\ Q M > / P d 

(Tvv)h = -g 2 m 2 Re[<0 M \T vv <f) m <f)l\0 M >]/P h (46) 

These two quantities are interpreted as the energy emitted when the atom is found in 
the ground state (disexcitation d) or in the excited state at t = +00 knowing that the 
atom was prepared in the excited state at t = — 00. 

When equilibrium is reached, the state of the atom is a thermal superposition. The 
probabilities of finding it in the excited or ground state are P + and P_ given in eq. (|38|) 
with P + + P_ = 1. Since Minkowski vacuum is a thermal distribution of Rindler quanta, 
one can approximate the state at equilibrium by 

\Aherm. > = P- > + P+ > (47) 

This neglects the dressing of the states due to high orders in H int but gives correctly the 
properties of the fluxes. The conclusions obtained using this naive state will be proven 
to be true to all orders in g 2 in Section |2^7. 



In the state eq. ( f47[ ) the energy flux is given given by the weighted sum of (Tw)f- 
and (T vv }^ + : 

(T VV (V)) therm. = P~(Tw)ip- + P+(Tw)ip+ (48) 

(This stems from the fact that the energy momentum operator changes the photon 
number by an even number and that the interaction hamiltonian changes the photon 
number by an odd number while changing the state of the atom). Hence all the matrix 
elements of Tyy we shall want to calculate can be expressed in terms of (Tyy)j where i 
stands for e, g, d and h. 

At this point, we remark that each each of these matrix elements (Tyv)i 

are acausal, 

for instance they are non vanishing in the left Rindler quadrant V < 0, U > 0. This will 



be discussed and interpreted in Section O . However the mean energies (Tvv)ii>- (where 



j = +, — , therm.) are causal. This follows from the following very general argument j4|. 
If V, U is separated from the trajectory of the atom U a (r), V a {r) by a space like distance 
then TyviYi U) commutes with Hi nt and one can rewrite the mean value as 

(T VV (V, U))^ = < iPjle+tf^^TwiV, U)e~^ dtH ^ |^ > 

= < ipjlTwiV, u)e +i $ dtHint e- i S dtHint \ipj > 

= <^\T vv (V,U)\^>=0 (49) 

Furthermore, by virtue of the conservation eq. (|39D , the zone where {TyviV, U))^. van- 
ishes can be extended to regions which are not space likely separated from the trajectory 
but to which V modes are unable to propagate. We note that a similar causality argu- 
ment applies in regions where (T vv (V,U))^ j ^ to guarantee that it only depends on 
Hint{r) for r's such that V(r) < V, i.e. that it only depends on the form of H int {r) in 
the past light cone of (V, U). 
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We also point out here a global property of the matrix elements (Tyy) e g which will 
have important consequences in the sequel. Namely that the Minkowski energy carried 
by (Tw)e is strictly positive 

/+oo 
dV{T vv ) e = g 2 m 2 < M \4>lH M <j) m \0M >/P e > (50) 
-oo 

since it is the expectation value of Hm (defined in eq. (|8|)) in a state which is not 
Minkowski vacuum. On the other hand the Minkowski energy carried by (Tyy) g vanishes 
identically 

/+oo r , 

dV(T vv ) g = -g 2 m 2 Re < Q M \H M (t>l<t>m^M > /P g = (51) 
-oo L J 

since H M \0 M > = 0. The same results are also true for (T vv ) d and {T vv ) h . 

In preparation for the next sections, we remark that all the matrix elements (Tyv)i 
(i = e, g, d, h) can all be expressed in terms of the following two functions (once the V 
term is dropped) 



C + {V) = < (WIW^JOm > = jdTG + (V,V a (T))e~ imT f(r) 
C-(V) = < O M |0OO0m|OM > = JdTG + (V,V a (T))e +mT f*(r) (52) 
where G+(V, V) is given in eq. (H). Indeed using eqs. (f|3|) and (filf ), one has 

(T vv ) e = 2{^){dvC-){dvCl) 
(T vv ) d = 2(^-) (d v C + ) (d v Cl) 

(Tvv) 9 = 2( 9 -—)Re[(d v C-)(d v C + )] = (-±)(T vv ) h (53) 

For these matrix elements of Tyv n ot to be singular the functions dvC + (V) and 
dyC-iV) must be regular. The function dvC + {V) can be expressed as 

d v C + {V) = ~l dr- 1 ar , f{r)e-^ (54) 

4tt J V — a L e aT — ie 

It can be singular only for V = where it takes the form 

dvC+( V) = -Lj iT _±_ f(T)e -^r ^JLJ dTe ^ f{T)e ^ (55) 

The last integral is finite if and only if /(r) decreases for r — > — oo quicker than e aT . 
Similarly if we had considered right movers, the condition for finiteness on the future 
horizon would have been sufficient rapid decrease of / for r — > +oo. Putting all together 
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the condition to not have singularities on the horizons is that /(r) decreases faster than 
e -o|r|_ This can be rewritten as 

Jdr^\f(r)\= jdt\f(r(t))\<oc (56) 

ie. the interaction of the atom with the field must last a finite Minkowski time. 

In order to obtain explicit expressions for the functions C± we reexpress f(r)e~ tmT 
in terms of its Fourier transform 

f{r)e- imT = J d\^-e- lXr (57) 

The normalization is 

J dr|/(r)| 2 = J d))—^— = T = total time of interaction (58) 



The regularity condition eq. ( p6|) is equivalent to having c\ be an analytic function 
of A in the strip —a < ImA < a. Hence we will not consider Lee models since they have 
c\ = for A < and thus singularities on the horizons. But in order that the behaviour 
of the uniformly accelerated two level atom be physically unambiguous, it is necessary 
that c\ be peaked around +m (the contribution of the negative frequency components 
of c\ should be negligible) and the golden rule probability of transition eq. fl34|) be 
recovered. For this to be the case T must satisfy T » rrT 1 and T » a" 1 . 

The first of these conditions is that f(r) be spread over a distance at least equal to the 
inverse frequency rrT 1 (the time-energy uncertainty condition). The second condition, 
which corresponds to T being greater than the euclidean tunneling time 2ira~ 1 |30f . 



is required for the probability P e to be linear in time and proportional to the Bose 
distribution N m . Both these conditions arise in the standard textbook calculation for 
the probability of transition in a thermal bath 

^ -99/", / x sin 2 (m — uj)T/2 

P e = 4g 2 m 2 / dwp{u) } ^— 59 

J [m — uj) z 

where p is the density of states (available photon states for disexcitation or the density 
of photons present for excitation). It is legitimate to replace [sin[(m — u)T/2}/(m — u>)] 2 
by 7rT5(m — id)/ 2 in the integrand provided T~ 1 dln p/du « 1. For disexcitation 
in vacuum this yields T » m~ l . In a thermal bath at inverse temperature /3, the 
Boltzmann distribution yields the additional condition T >> (3, hence T >> a -1 in the 
uniformly accelerated case. 

We now express the operator 4> m and the functions C± in terms of c\. 

00 ^ a\,R * f° ^ a U,R * 
dX , =c A + / dX =c x 

aMttIAI aMttIAI 



^/47rA(e 7rA /a - e -TV«) 
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+00 



d\c\ 



+00 



d\c\- 



C-{V) 



-00 
+00 



d\cX 



y / 47rA(e 7rA /« - e -^/a^ 


1 

4vrA 


\n x + l)(aV)- tX / 




1 

1 



,7rA/2a 



47rA(e 7rA / 



+00 I 

d\cX- — - 
A 47rA 



wX/a \ 

n x {aV) iX/a 6(V) + n x e wX/a \aV\ iX/a 6(-V) 



(60) 



where we have used eq. (|T5|) for the expression of Lp XjM (V). And n x = l/(e 27TX ^ a — 1) is 
equal to (see eq. (|i~9"D) 



"A 



N x =P 2 X for A > 



-{N w + 1) = -af X] forA<0. 
The probability P e to excite can also be written in terms of c x 



2 2 



+00 



dX 



\cx\\ 



4vrA 



^A 



(61) 



(62) 



As one picture is worth a thousand words we take a particular form for c x such that 
all the integrals above are gaussian and can be evaluated explicitly 



A 



CA = jD _ e -(A-) 2 TV2 (1 _ e 
m 



-2-wX/a 



(63) 



where D is a normalization constant taken such as to verify eq. (|58|). 

We shall give throughout the text the exact expressions followed by the approximate 
expressions which are valid when T satisfies the condition discussed above T » rrT 1 
and T >> cT 1 as these last are physically relevant and are particularly easy to read and 
understand. The approximate expressions are preceded by the symbol ~. We conclude 
this section by the value of the switch off function / (see figure 2) 



f(r) 



7T 



D c -r 2 /2T 2 
-1/4 -r 2 /2T 2 , 



I T \ _ e -27rm/a e j27rr/aT 2 e 27r 2 /a 2 T 2 ^ 1 _ ^ ' 

mT 2 
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mT 2 amT 2 ' 



l + N m (l-e l2 * r/aT2 )) 



(64) 



where the constant D takes the form D ~ 2 1 / 2 7r 1 / 4 T(iV m + 1). Eq. (|4|) shows the 
almost gaussian character of the switch off function whose width is T. The plateau of 
the gaussian gives a good approximation of the steady state regime which we which to 
study. 



2.5 Fluxes and Particles to Order g 2 During Thermalisation 

We briefly sketch the main results of this section. During thermalization a steady flux 
of negative Rindler energy is emitted: (T vv ), ~ —g 2 m 2 N m /2. This is understood 
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from the isomorphism || with the thermal bath: as the atom gets exited it absorbs 
energy from the thermal bath, thus the minus sign. The transcription of this flux to 
Minkowski quanta is more subtle. Oscillatory tails in the Rindler flux are enhanced by 
the jacobian that converts from Rindler to Minkowski energy with the net result that 
positive Minkowski energy is emitted. In the Minkowski description the origin of the 
steady negative flux is due to a "repolarization" of the atom corresponding to the fact 
that the probability of finding the atom in its exited level decreases with time. This 
repolarization is similar (CPT conjugate) with that which occurs when negative energy 
is absorbed by an inertial detector . 

We shall discuss both the adiabatic switch on and off presented in the previous section 
(to reveal the oscillatory tails) and a sudden switch on and off (to display the properties 
in the stationary regime). 

We start with the adiabatic switch on and off. In terms of the function c x introduced 
in eq. ( j63| ) the mean energy radiated by the two level atom initially in its ground state 
is 

(Tw(v%_ = -9 2 m 2 J d\j 'dX'c x cl,j^(n x + h xl )e-^ x - x > 

" iV m ^ 7 ^[(iV m + l)cos(27rt;/aT 2 )-iV m ] (65) 



2 m Tr 1 ^ 
As announced it carries negative Rindler energy: 



-hx g 2 m 2 r+oo 



dv{T vv (I + )),_ = — / d\\c x \ n x 

~ --g 2 m 2 N m T = -mP e (66) 

which is equal to the probability to be found excited times the absorbed Rindler energy 
— m. 

The total Minkowski energy radiated is 

r+oo r+oo 

(Hm) c = I dV(T vv (V))^_ = 1^ dve™(T vv (v))^_ 

~ +^g 2 m 2 N m Te aT °(l + 2N m ) = +mP e e aT °(l + 2N m ) (67) 

where e aT ° is the mean Doppler effect associated with the window function /(r). We 
define it by 

J dve~ av e- v2/T2 cos(2v™/aT 2 ) = -e^T/vr 1 / 2 (68) 

The Minkowski energy is positive (as it should be), whereas the Rindler energy is neg- 
ative. The flip in sign is due to the effect of the transients around v = aT 2 . Indeed 
whereas the transients are negligible upon computing the Rindler energy eq. (0), upon 
computing the Minkowski energy they are enhanced by the jacobian dv/dV and give 
rise to the flip in sign. (Note that this sign flip of the Minkowski energy versus the 
Rindler energy can be conceived as arising from the imaginary part of the saddle point 
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of eq. d68|): v sp = —aT 2 /4 + i7r/a and is therefore on the same footing as that the flip of 
frequency which leads to a non vanishing (3 coefficient, see ||30||). The additional factor 
1 + 2N m in eq. (|6T| ) and the difference with eq. (0) by a factor Ta comes from the 
inherent ambiguity in defining e aT ° as the mean Doppler shift associated to the switch 
function /(r). 

We insist on the fact that the total Minkowski energy radiated can also be expressed 

as 



/+oo 
dV(T vv ) 
-oo 



(69) 



because of eq. (|51|). So the Minkowski energy can be conceived as coming from (Tyy) e 
only. When expressed in this fashion the integrand in eq. (|69|) is strictly positive but 
located essentially in the region V < (this is shown in Part 3). (We remark that 
eq. (§) can be viewed as a rewriting of eq. (|69|) but with the integration taken in the 
quadrant V < 0). The (Tw)g term restores causality and localizes all the energy in the 
transients. 

Another case of interest is the golden rule limit for which c\ = 2tt5(X — m) corre- 
sponding to /(r) = 1 for all r. In this case their is a constant negative flux for all 
V > 0. The transients are located on the past horizon V = where they consist of a 
singular positive flux j7|. Rather than this case we now analyse the case where the time 
dependent coupling is /(r) = 6{t)6{T — r). With this time dependence the transients 
are singular and will not be studied (this divergent behavior is already present for an 
inertial detector with the same switch function and has nothing to do with the presence 
of a horizon). On the contrary, the steady part is easily computed and corresponds ex- 
actly to the intermediate values (— aT 2 « t « aT 2 ) found in the adiabatic situation 
described above in eq. (p5|). The reason for which we shall now belabour this case is 
that it shows explicitly the relations between the flux emitted and the transition rate 
(not only the relation between the probability and the total Rindler energy as in eq. 

©)■ 

The probability of spontaneous emission is given by 

r T f-T 

Pe{T) 



2 m 2 / dn / rfr 2 e~ im ( T2 - T1 )(0(r 2 )0(r 1 )) 
Jo Jo 

2„ 



g mN m T 



(70) 



The second line contains the golden rule result valid when aT — > oo with g 2 T finite. It 
is useful to introduce the rate of transition, the derivative of P e {T): 

fT 



p = dP^T)_ = 2 m 2 2Re 

e\ ) dT 

~ \g 2 mN m 



dre- im(T - r) (0(T)0(r)) 



(71) 



This rate is related to the (steady part of) the stress energy tensor. Indeed one finds 



g 2 m 2 2Re 



T 



drie 



([0(t 2 ),T to (T)]_0(t 1 )) 
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g 2 m 2 2Re 



dTe- im{T - T \id v <P(T)<p(T)) 



-mP e (T) + g 2 m 2 2Re 



te 



-imT 



<0(T)0(O)> 



(72) 



The first equality follows straightforwardly from the expansion of the evolution oper- 
ator z~ % l HintdT in g 2 . The second equality is obtained using the commutator relation: 
[4>{t2), T w (ti)]_ = id v (f)5{ri — r 2 ). The third equality follows by integration by parts. The 
final result contains a steady part proportional to —mP e (T) which tends to — ^g 2 m 2 N m 
in the golden rule limit and an oscillatory term (which is exponentially damped if a 
slight mass is given to 0). The steady piece simply indicates that to an increase of the 
probability to make a transition corresponds the absorption of the necessary Rindler 
energy to provoke this increase. 

We now turn to the Minkowski description of this steady piece. We first rewrite these 
expressions in terms of the Minkowski basis e~ l /v4ttu> (see eq. (^)). The probability 
of transition eq. (f70|) reads 



Pe{T) 



1 1 

g m 



du 



dr e 



y/Anu 



POO 

/ duP e ^(T) 
JO 



(73) 



Similarly the transition rate eq. (|7T|) becomes 



P e (T) = g 2 m 2 / du 2Re / dTe~ im{T - T) 
Jo Jo 



±(e aT -e aT ) 



o 



du Pe,u>(T) 



And the total Minkowski energy is given by 

r+oo 



(Hm(T)) ( 



dve° 



POO 

/ du uP e>ul (T) 
Jo 



(74) 



(75) 



(Where in the first equality the integral is only over region of positive V since by causality 
the mean energy is unaffected in the other quadrant, see eq. (|4*9"|)). The second equality 
follows from the diagonal character of the energy operator Hm see eqs. (|50|), (|5T|). The 
positivity of {Hm{T)) £ is manifest since all the P e ^{T) are positive definite. Nevertheless 
the time derivative of (H(T)) e is negative, within the steady regime, 



d(H M (T)) t 
dT 



du uP e ^(T) 



^ T \T vv (v(T))\ 



-me 



av(T) 



P e {T) + oscillatory "damped" term 



(76) 



d{HM)/dT negative implies thus that, for large u (since P e (T) > 0), some P eyU) are 
negative. This corresponds to a "repolarization" since all the P e ^ are positive definite 
and vanish for r < 0. This repolarization is exactly the inverse process of the absorption 
of negative energy by an atom described in |34]] . 
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2.6 Fluxes and Particles to Order g 2 at Equilibrium 

Before studying the equilibrium situation it behoves us first to consider the flux emitted 
by an atom that makes a transition from excited to ground state. 
The mean energy emitted when the initial state is > is 



(Tw}^ + = g 2 m 2 Jd\ Jd\'c x cl,—^(n x + nx> + 2)e 



i(X-X)v 

(47 



~ ^(iV m + l)e-" 2 / T2 [l-iV m {cos(27ri;/aT 2 )-l}] (77) 



2^ 

and the total Rindler energy radiated is 

g*m 
1 



2^2 



dv(T vv (v)}^ + = I d\\cx\ 2 (n x + 1) 



~ -g 2 m 2 (N m + l)T = mP d (78) 

In the example for which the time dependent coupling is f(r) = 6(r)9(T — r), the 
relation between the derivative of the probability P d (T) and the flux (T TO )^ is 

(T VV (T))^ = +mP d (T) + oscillatory "damped" term (79) 

The sign in front of Pd{T) is now positive (contrary to the one in eq. (|72|)). Disexcitation 
consists in emitting the energy stored in the atom. 
The total Minkowski energy emitted is 

p+co q 2 vt? 

J o dV (T vv )^ + ~ —{Nm + l)Te aro (2N m + 1) = mP d e aT °(2N m + 1) (80) 

For the disexcitation, the integrated Rindler and Minkowski energies have the same sign 
and are related by the mean Doppler shift e ar ° times (2N m + 1). 

We now turn to the thermal equilibrium situation. We recall that the energy radiated 
is the sum of the fluxes emitted when the atom is initially in its ground state and when 
the atom is initially in its exited state weighted by their initial probabilities. Hence one 
has 

(Tw/ 'therm. = P~ (Tvv) + P+ (T vv ) ^ + 

~ -mP_P e + mP + P d = (81) 

The steady fluxes cancel exactly each other because at thermal equilibrium P± satisfy 
eq. (P8"|). This is Grove theorem in g 2 Only the oscillatory transients remain. 

They read 



(T m ) therm . = ^^--JdxJdX'^ 



(4tt) 2 

[N m (nx + n x , + 2) - (N m + l)(n A + fix)] e~^- A > 

2 2 

9 m :N m {N m + l)e~ v2/T2 \l - cos(2t™ /aT 2 )] (82) 



Arc 
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To illustrate the positive transients, we have plotted (T vv ) therm in figure 3. 
The total Rindler energy emitted is 

dv {T vv ) therm . = ^ 2Nm + 1 J d Mcx\ 2 (N m -h x ) 

~ ^-N m (N m + l)^- (83) 

It tends to zero as the time of interaction T tends to oo i.e. as c\ tends to a 5 function. 
(In this case, the two level atom tends a Lee model. This can be seen in eq. ( p3|) where 
the negative frequencies are exponentially suppressed.) 

However, the total Minkowski energy increases with the interaction time T and is 
given by 

r+oo r+oo r+oo 

dV(T vv ) therm . = P_ dV(T vv ) f _ + P + dV(T vv )^ + 

~ m(P_P e + P + P dtV )Te arQ (2N m + l) (84) 

The Minkowski energy of the two fluxes coincide, by virtue of eq. ( |38| ) and sum up. This 
result is what one might have "naively" guessed : The total energy is the integral over 
the interacting period of the rate of transition times the varying Doppler shift times the 
energy gap m. 

We now go to all order in g to prove that this emission of Minkowski quanta is not 
an artefact of the second order perturbation theory. The forthcoming section can be 
skipped by the reader mainly interested by the study of vacuum fluctuations and the 
black hole problem. He can go to Part 3 directly. 

2.7 Fluxes and Particles to All Order in g 

We use the exactly solvable model (RSG), used by Raine, Sciama and Grove @-||, 
to prove that one does recover, to all order in g, that every quantum jump of the 
accelerated oscillator, in thermal equilibrium in Minkowski vacuum, leads to the emission 
of a Minkowski quantum. Hence the rate of production of the Minkowski quanta is simply 
the rate of internal transitions of the oscillator. But, as in second order perturbation 
theory, these quanta interfere and their energy content is found at the edges of the 
interacting period only. This is due to the complete neglection of the recoils of the 
oscillator. (Upon taking into account the recoils by giving the oscillator a finite mass, 
i.e. by quantizing the position of its center of mass, one proves that the Minkowski 
quanta no longer interfere after a short time (a few l/o)[32|). We conclude this section 



by giving a model independent proof that the stationary thermal Rindler equilibrium 
corresponds to a production of Minkowski quanta. 

We first recall the main properties of the RSG model and then analyse the particle 
content of the emitted fluxes. 

This system consists of a massless field coupled to a harmonic oscillator maintained 
in constant acceleration. Its action is 

S = Idtd.r - \{d t (j)f - (d x (f)) 
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+ dr 



(d T q) - m 2 q 2 + e(d T q)(j) 



(85) 



where X M (r) is the accelerated trajectory eq. (p6|) and e = g\ / 2m is a rescaled coupling 
constant. Since this action is quadratic, the Heisenberg equations are identical to the 
classical Euler Lagrange ones. They read: 



d u d v <f> = ~9(V)5(p-l/a)d T q 



(86) 



d 2 q + m 2 q = -ed T (j){X^{r)) (87) 

The left part of the field (i.e. for V < 0) is, by causality, identically free. And, for 
V > 0, on the left of the accelerated oscillator trajectory, the w-part of the field only is 
scattered. There the general solution is 



[u,v] 



(u) + (j){v) + -q(v) 



+oo 



q( v ) = q ( v ) +i I d\ ip x e iXv [<f>xjt, v + <j>\,R, u ] 



(88) 



(89) 



where (p{u) and 0(f) are the homogeneous free solutions of eq. (|36"D; where the /operator 
4>x t R )V is defined by 



1 [ 6(\)a KR + 6{-\)al KR 



>\,R,v 



4vr|A| 

^a similar equation defines 0A,i?,,«); where ip\ is given by 

eA 



(90) 



m? — A 2 — ie 2 X/2 



(91) 



and where q(v ) is a solution of 



d 2 q + m 2 q H d T q = 



(92) 



The two independent solutions of eq. (|92|) are exponentially damped as r increases. 
Being interested by the properties at equilibrium, we drop q(v ) from now on. Then, the 
remaining part of q(v) is a function of the free field only. Hence, in Fourier transform, 
eq. (|88D reads 



6 6 



(93) 



The second term in eq. (|9~3"D mixes m and v modes. It encodes the static Rindler polar- 
ization cloud (see ||) which accompanies the oscillator and carries neither Minkowski 
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nor Rindler energy. In order to simplify the following equations, we drop it and mul- 
tiply the other scattered term by two for unitary reason -see below. (By a simple and 
tedious algebra, one can explicitly verify that this modification does not affect the main 
properties of the emitted fluxes). Then eq. (p3|) becomes 

4>\,r,v = 4>\,R,v(l + iei)\) (94) 

It is useful, for future discussions, to introduce explicitly the scattered operators dx,R, 
and the scattered modes <p\ t R(v) 

d\,R =< <P\,r\4> >= a\ )R (l + ieipx) (95) 



fl A„R> • 



1 Px,r{ v ) = 

whereupon the scattered field operator <fi(v) may be written as 

~ roc 

<j){v) = dX [dx,Rf\,R + h.c] 
Jo 



1 + ietJ) X )<p X)R (v) (96) 



poo 

/ dX [a x , R (px,R + h.c] (97) 
Jo 



It is now straitforward to obtain the scattered Green function and its Rindler energy 
content. If the initial (Heisenberg) state is Minkowski vacuum the w-part of the scattered 
Green function is, for V, V > 0, 

G+{v,v') = <O m |0WK)|Om> 

dX\l + ieipxl 2 {Pa<P\,r(v)(P\,r( v ') + a lvxA v )^\A v ')) 
= G+(v,v') (98) 

where G + (v, v') is the unperturbed Minkowski Green function and where we have availed 
ourselves of the identity (see eq. (|9TD) 

\l + ie^ x \ 2 = l (99) 

This unitary relation expresses the conservation of the number of Rindler particles. 
Indeed there is no mixing of positive and negative frequencies in eq. (|95"D ; in other 
words, the /3-term of the "Bogoljubov" transformation eq. ( p5| ) vanishes. 

The identity of the Green functions in eq. (^) proves that, once the the steady 
regime is established, no flux is, in the mean, emitted. This is Grove theorem || 0. 

We now examine how this stationary scattering of Rindler modes is perceived in 
Minkowski terms. The Minkowski scattered modes <fx,M are given by 



</?a,m(1 + iea\il)x) - ieaxM\¥*-x,M 



23 



<P-\,M 



V-a,m(1 - ieplip-x) - iea\(3\ip-\(p XtM 

= a-\ip-X,M + P-X<P*\,M ( 101 ) 

where < A < oo and where we have introduced the scattered Bogoljubov coefficients: 



a x = 


1 + ieotxipx 


h = 


-ieax(3\ip*x 


&-X = 


1 + iePl^x 


0-X = 


-itaxMx 



(102) 



One verifies that the unitary relation is satisfied: \&x\ ~ \Px\ = 1- The fact that 
the (5 are different from zero indicates that each couple of jumps of the oscillator (the 
absorption and subsequent emission of a Rindler quantum) leads, in Minkowski vacuum, 
to the production of two Minkowski quanta. The member (p-x,M is emitted when the 
oscillator absorbs a rindleron and jumps into a higher level and the other one, (px,M is 
emitted during the inverse process. This is manifest in the mean energy flux: 



(Tvv) 



iw^dvdv.il^vwv') - <p(y)(fi(y')\) 



+oo 



dX \px\ 2 \d V ¥x,M\ 2 + Re ®xf3*xd V ¥x,Mdv¥-x, 



..M 



(103) 



whereupon the total Minkowski energy is 

r+CO 



(H M ) 



dV(T vv ) 

oo 

+°° ~ 2 ~ 2 f+°°dV 



dx x(\p x \'+\p-xn 



■3C 



2vr a?\V + ie\ 2 



(104) 



since the integral of the second term vanishes. 

Exactly as in second order perturbation theory, there is a steady regime during which 
all the emitted quanta interfere destructively leaving no contribution to the mean flux 
(see eq. (^)). But all non diagonal matrix elements will be sensitive to the created 
pairs. This is also the case for the the total energy eq. (|104|) since being diagonal in u it 
ignores the destructive interferences (the second term of eq. ( |103| ) whose role is to make 
the mean flux vanishing during the steady regime). 

In order to prove that eq. (|104|) corresponds to a steady production of Minkowski 
quanta during the whole interacting period At = T (infinite in eq. ( |104| )) we evaluate 
how many quanta are produced. (Contrary to the energy, the total number of Minkowski 
quanta is a scalar under the Lorentz group, hence not affected by the exponentially 
growing Doppler shift present in the energy) 



r+oo 

(iV(Ar)) = / du<0 M \ala u \0 M > 
Jo 

r+oo 

= / dui< M \ala w \0 M > 
Jo 



+co r+oo 



duj / dX \ 1XiUJ (Ar)\ 2 \i3x\ 



(105) 
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where \0m > is the scattered (Schrddinger) statefj The are related to the a XM by 
(see eq. (fjjD) 

/•OO 

a w = / rf^ 7A,u;(Ar)a Ai Af (106) 



where jx,oj (At) takes into account the time dependence of the coupling. As shown in 



30||31| 7^,0; (At) is non vanishing only for the u which enter into resonance with the 



oscillator frequency m during the interaction period Tj < r < tj = Tj + T. When these 
frequencies belong to 

w< = me"" T ' < u < me' aT f = u f (107) 
7A,u;(Ar) may be replaced by 7a,u; (given in eq. flTBD). Hence iV(Ar) reads 

(N(Ar)) = f"£-r<*M 

A T r+oo 

d\ \Px\ 2 (108) 



27T J-oo 

The total energy emitted obtained from eq. (|108|) is 

J / da; 



(^m(At)) = /'_/ dA|/5 A | 



^e~ aT m / dA |/3 A | 

VfdV 1 



OO 

' OO 



m / dA |/? A | (109) 



Vj 2vra 2 V 2 

in perfect agreement with eq. ( |104| ) if the frequency width of the oscillator in small 
compared to m. The rate of production (eq. ( |108j ) divided by At) is (small width limit) 
e 2 ctmfim which is the rate of jumps for an inertial oscillator in a bath at temperature 
a/27r. Therefore the number of Minkowski quanta produced by the thermalized oscillator 
equals the number of internal jumps. 

We now generalize these results to an arbitrary linear coupling. We believe that it 
can be generalized, using the same type of argumentation, to nonlinear couplings as 
well. The proof goes as follow. Any scattering of Rindler quanta by an accelerated 
system which leads to a thermal equilibrium during a time much larger than 1 /a can be 
described as in eq. ( p5|) by 

a\, R = S X \>a x >,R (110) 

where repeated indices are summed (or integrated) over and where the summation over 
A' includes both u and -u-modes (as in eq. (p3|)). The matrix S satisfy the unitary 
relation 

S XX "Sy,y = 5 XX > (HI) 

5 The simplest way to obtain this state is to find the scattering operator U such that &a,m = U^cl X: mU 
where a x . M = (<f X ,M\4>)- Tnen \®m >= U\0 M >■ 
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which express the conservation of the number of Rindler quanta since S\ t \' mixes positive 
Rindler frequencies only. It is convenient to introduce the matrix T (from now on we do 
not write the indices) 

S = 1 + iT 



i a,_r; a \,L 



which satisfies 

2ImT = TT+ 

We introduce also the vector operator b = \ a\ t R] ax,L] • 
written as 

~b = Sb 

where S has the following block structure 

/ 1 + iT \ 

10 

1-fft 

V 1 / 



(112) 
(113) 

Then eq. ( |110|) can be 
(114) 



(115) 



since the u and f-modes on the left quadrant are still free. On the other hand, the 
Bogoljubov transformation eq. (|T8|) reads in this notation 

c = Bb (116) 
where c = (a\ t M', o>-\,m'-> a A,A/> a -A,Af) an d where B is 



B 



fa -0 \ 

a -p 

-0 a 

V -p a J 



(117) 



the diagonal matrices (in A) a and P being taken real. The scattered Minkowski operators 
are given by 



c = BSB~ l c ={S + B 



S,B 



-i 



c = S M c 



(118) 



Since S and B do not commute, Sm has non diagonal elements which encode the pro- 
duction: 



S 



M 



/ 


Oil 








-A 























A 









V 


-fh 








a! 2 


) 



(119) 



where the a j3 are given in terms of T by (see eq. (|102[) ) 

di\ = 1 + iaTa 



Pi 
5 2 

a 2 



—iaTP 
1 + i/3T f /3 
iPTa 



(120) 



QED 



26 



3 The Conditional Values of the Energy Momentum 
Tensor 



3.1 Introduction 

In Part 2 we analysed the mean energy radiated by the atom and showed, in Section 
2.4, how it can be decomposed into two contributions (Tyy(U, V)) e and (Tw(U, V)) g 
according to the final state of the two level atom. These correspond to the energy emitted 
if the atom is found in its excited or ground state at t = +oo. This decomposition is 
valid for all points U > U a (V) where U a (V) is the trajectory of the atom, i.e. in the 
future of the atom since for our massless field energy flows along V = constant. 

In this Part, we generalize this decomposition of the energy density into energy densi- 
ties correlated to the final state of the atom for all U, V. When U < U a (V) this describes 
the energy momentum of the field configurations which will give rise to excitations of 
the atom. In order to explicitize this proposition, we proceed as follow. First, the gen- 
eralization is introduced formally by introducing projectors that specify the state of the 
atom at t = +oo. We shall then see that this construction gives rise to nondiagonal 
matrix elements of Tyy which are complex. We need therefore to discuss their physi- 
cal meaning. To this end one introduces an additional quantum system coupled to the 
operator Tyy. The picture that emerges is then clear: to first order in the coupling, 
the modification of the wave function of the additional system is governed by these non 
diagonal matrix elements. This is presented in summary fashion in Section 3.2 and in 
more detail in Part 5 which is entirely devoted to a general discussion of the procedure 
leading to these conditional values of operators. In the present Part, we discuss mainly 
the properties of these conditional energy densities. 

In section 3.3, the conditional values of the energy distribution correlated to the 
transitions of the accelerated atom are described and interpreted. In Section 3.4 we 
discuss a generalization of these conditional values which does no longer refer to the 
transitions of the atom and which finds important application in the black hole problem. 



3.2 The Conditional Energy Correlated to a Transition of the 
Accelerated Atom 

In Section 2.4, we had rewritten the mean energy emitted on the left of the accelerated 
trajectory (i.e. U > U a (V)) as 

(Tyy)^_ = P e (Tyy) e + P g (Tyy) g (121) 



using the expression, eq. ([31), for the > at t = oo and the probability of transition 
P e given in eq. fl34). In order to generalize this decomposition, we first rewrite it in 
Heisenberg representation by introducing the projectors Il + = |+ >< +| and EL = 
|— >< — | onto the excited and ground state of the atom. 

In Heisenberg representation, the state of the system is > = \0m >\— > and the 
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projector is a time dependent operator given by 

n + (t) = e i J-°°* aiat Il + e~ i J-°°' kHint (122) 
whereupon the probability to be found in the excited state at t — +00 is written as 

P e = (^_|IT+(t = +oo)|^> (123) 
Similarly, the probability to be found in the ground state at t — +00 is 

P g = (V-iII_(t = +oo)|V>-) (124) 

The conservation of probability P e + P g = 1 is realized through the completeness of the 
projectors 11+ (t) + IL_(t) = /. 

The conditional energies can be now obtained by decomposing the mean using the 
projectors 11+ (t) 

{t vv {u,v))^_ = (V-|[n+(+oo) + n_(+oo)]T vv (c/,y)|^-) 

(^\U + (+oo)T vv (U,V)\^) (^\U_(+^)T VV (U,V)\^_) 

(^_|n+(+oo)|^_) 9 (^_|n_(+oo)|^_) 

= P e (Tw(U,V)) e + P g (T vv (U,V)) g (125) 

When U > U a (V) , the matrix elements (Tyy(U, V)) e and (Tyy{U, V)) g are the expres- 
sions obtained less formally in eq. (f|^). When U < U a (V), these matrix elements are 
the desired expressions of the energy density if the atom shall be found at t — +00 in 
the excited (ground) state. 

The normalization in eq. ( 125 ) is chosen so that the mean value is expressed as the 
probability of making a transition times the conditional value exactly like in the usual 
conditional probabilities. We shall see in the sequel that this decomposition into proba- 
bility to end up in the excited or ground state times the conditional value automatically 
occurs in physical processes. 

Two important properties of the conditional values for U < U a (V) should be noted. 
First 

(Tyy(U < U a (V), V)) e = -^-(Tyy(U < U a (V), V)) g (126) 

since (T VV (U, V))^_ vanishes identically for U < U a (V) because the interaction with 
the accelerated atom has not yet perturbed Minkowski vacuum. Secondly, (T vv ) e is 
complex. This can be seen from the explicit expression 

(Tw(U<U a (y),V)) e = ^(tfj.\e l I dtH ^U + e- l I dtH ^T vv (U,V)\tlj.} 

= ^(0M|^mTvy(c/,i/)i0M}^c;(y)c*(y)(i27) 

where we have used eq. ([IT]) and eq. fl5"2"|). 
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Both the real and imaginary part of (Tyy) e have physical meaning and intervene 
in physical processes. To prove this fact one should introduce an additional quantum 
system because these matrix elements have meaning only in quantum mechanics. 

For definiteness, we take the additional system to be a quantum oscillator sitting at 
X = Xq and coupled to Tyy by the interaction hamiltonian 

J dtH osc . = J dtg vv (t)p(t)Tyy(t,x ) (128) 

where p(t) is the momentum conjugate to the position q(t) of the oscillator and g (t) 
is a switch function with the correct Lorentz variance. The initial state of the oscillator 
is \osc. >. The state of the entire system (i.e. field + two level atom + oscillator) is 
thus |^ >_ = \ip >_ \osc. >. 

Then to first order in g vv , in the interacting picture, the mean position of the 
oscillator at t = oo is given by 

(<?(£ = oo))^_ = < osc.\q(t = oo)\osc. > 

+ J dtg vv {t)< osc.\ —i[q(t = +00), p(t)]-\osc. >{T vv (t, x ))</,_ 

(129) 

That is, the mean change of the position is driven by the mean value of Tyy(t, Xq) in 
the state >_. It corresponds to the classical response of q{t) to a driving force. 

But, one can also investigate the correlations among the oscillator state and the atom 
by can asking more detailed questions such that: what is the "mean" (for the use of 



this word see discussion before eq. (|4)) position of the oscillator when the two level 
atom is found in its excited state? The answer is the conditional value of q obtained by 
decomposing the mean according to the final state of the atom at t = 00 

(q(t = +oo)>*_ = P e {q(t = +oo)) e + P a (q(t = +oo)) 9 (130) 

where the conditional value (q(t = +00)) e is given by 

""■^ ^lat?" <i3i) 

To first order in g vv , this conditional position of the oscillator is 

(q(t = +oo)) e = < osc. \q{t = — oo)\osc. > 

dtg vv (t) (osc. \ —i[q(t = +00), p(t)]-\osc.) Re[(T vv (t, x )) e } 

+ / dtg vv (t)(osc.\-i{q(t = +00), p{t)} \osc.) lm[(T V y(t,x )) e ] 



(132) 

Hence both the real and imaginary part of the conditional value of Tyy control the 
modification of the mean conditional position. Note that Re{Tyy) e for U > U a (V) and 
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for U < U a (V) enter exactly in the same way in the integrals giving rise to (q(t = +oo)) e 
as the mean value (T vv (t, x ))^_ drove the mean q in eq. (|129|) . The imaginary part of 
(Tw)e appears in an unusual way through an anticommutator which depends explicitely 
on the state of the oscillator. In quantum mechanics therefore, by coupling an additional 
system to the operator Tyy, one can isolate, in a well defined manner, the energy content 
of the emitted particle correlated to a transition of the atom as well as the energy content 
of the vacuum fluctuations that shall induce the transition of the atom at later times. 

This procedure wherein an external quantum system is introduced to reveal the 
physical significance of matrix elements like (Tw)e,g will be displayed in more details in 
Part 5 and put in parallel with the treatment of Aharonov et al.[0- The same procedure 
will also be used in the black hole situation when evaluating the conditional value of the 
metric correlated to a particular final state of the radiation. 

3.3 The Properties of the Conditional Energy 

Having indicated by an example how both the real and imaginary parts of (Tyy) e inter- 
vene in physical processes we now display the properties of the conditional values. Since 
(Ty V (t,z))g = ({T vv )ip_ - P e (T vv ) e )/P g we shall discuss (Tyy) e only. 

In order to obtain exact expressions for this matrix elements, we use again the c\ 
introduced in eq. (p3"|). We give now the three expressions for {Tyy) e : three because 
one finds different expressions for V > 0, U < U a (V) and for V > 0, U > U a (V) (i.e. 
before or after the interaction occurs), and for V < all U's. 




(133) 




ml 



1 - 



iv + 27r I a 



mT 2 



2^T e 




(134) 





m(N m + l) iv L +7r/a l2 4 

1 — ■ P 



(135) 
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where the last two equalities in eqs. ( |133| , |134| , |135| ) represent the exact expressions if c\ 



is given by eq. (|63|) and the approximate expressions valid for T >> m , T >> a 
Cq is a constant equal to 



Co = (N m + 1) 



'I 71 ) c -2^m/a c ^/a^M 2lX 



amT 2 amT 2 ' 



1 (136) 



These fluxes are presented in figure 4. 

We now present the complementary Rindler and Minkowski properties of these con- 
ditional values of T vv . 

The Rindler description is that used by a uniformly accelerated observer in the same 
quadrant as the two level atom. It is best understood by making appeal to the isomor- 
phism of the state of the field in the right Rindler quadrant with an inertial thermal 
bath. 

By getting excited the two level atom has selected that the thermal bath contains at 
least one particle in the mode created by <$ m . Furthermore since energy flows along the 
lines v = est, (T VV (U < U a ,V > 0)) e is centered around v = with at spread Av = T. 
It carries a Rindler energy obtained by integrating eq. ( |133j ) 

dv{T m (U <U a ,V> 0)) e = /rfA ' C ;'|^- +1) - + 1) (137) 

The factor N m + 1 takes correctly into account the Bose statistics of the field since eq. 
(|137|) corresponds to evalutating (n 2 ) / (n) in a thermal distribution. 



Then, by getting excited the two level atom absorbs one quantum and the residual 
energy on the future of the accelerated trajectory U = U a (V) is (see eq. ( |134j )) 

/ dv(T„(U >U a ,V> 0)) e = /!^ |C ff - mN m (138) 

We now consider what is "seen" by a uniformly accelerated in the left Rindler quad- 
rant (i.e. what is the nature of the correlations between the transition of the atom and 
an additional system uniformly accelerated in the left Rindler quadrant). Before the 
strict correlations between the left and right quadrants, see eq. (p5|), one expects that 
to the (N m + 1) Rindler quanta on the right correspond (N m + 1) Rindler quanta on 
the left. This can be obtained formally by considering the Rindler energy operator Hr. 
Since Minkowski vacuum |0m > is annihilated by Hr (see eq. (|J) the Rindler energy 
in the left quadrant is equal to the energy in the right quadrant. Indeed integrating eq. 
(|135|) and using the relation fi\{n\ + 1) = n 2 x e 2nX ^ a yields 



J dv L (T VLVL (U <U a ,V< 0)) e = J dv(T vv (U <U a ,V> 0)) e (139) 

The symmetry between the left and the right Rindler quadrants results in (T VV (U, V < 0)) 
being centered around vl = with the same width Avl = T. Thus (Tyy(U,V < 0)) e 
carries also the Rindler energy of N m + 1 Rindler quanta and is almost exactly the sym- 
metric of (Tw(U < U a (V),V > 0)) e except for small transient oscillations present for 
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V > (see the explicit expressions eqs. ( |133| , |135| ) and figure 4). We also note that 
(Tyy(U, V < 0)) e is real whereas (T VV (U < U a (V),V > 0)) e is complex. This results 
from causality and can be proven in complete generality by making appeal to a reason- 
ing similar to that in eq. (f49|) . It will have important consequences in the black hole 
problem, in Section ^4.5| . 

The Minkowski description, i.e. that used by an inertial observer, is best understood 
by rewriting the conditional value of T vv in terms of the <^a,m(U) modes, eq. (|T5|), 

j (P'TT? ft I I 

(T vv (U<U a ,V)) e = a2y2 y pe jd\jd\'c* x cx—J\\'nx(n x + 1)^, M ^_ A ,, M 

1 m(N m + 1) i . TT 

1 + ZZZ^i H-aV - ze) - — — ) 



a 2 V 2 2^TC v maT 2 v ' maT 2 

x (1 l — \n(-aV - ie) ^ y-[M-aV-ie)f/a^ (14Q) 

maT 2 maT 2 

The ie defines ln(— aV — ie) as In \aV\ for V < and as In \aV\ — in for V > 0. Upon 
taking the limit e — > no singularity occurs. In fact (Tyy(U < U a ,V)) e given in eq. 
(|140|) vanishes for V = 0. This is an accident due to the particular form of c\ chosen 



in eq. fl63|) (it has zero's for A = ina, n = .., —1,0, 1, ..). But, from the expression for 
C±(V = 0) given in eq. fl55p , it results that the generic behaviour of Tyy is to stay 
finite as V — > 0. In more physical terms this corresponds to saying that the Minkowski 
vacuum fluctuation that induces the transition straddles the horizon with no clear cut 
separation between the pieces in the left and right quadrants. 

Notice how the ie prescription which encodes the analyticity of the modes (f\ t M m 
the lower half complex plane now encodes the vanishing of the integral 

/+oo 
dV (Tyy(U<U a ,V)) e = (141) 
-oo 

by contour integration. The vanishing of this integral can also be seen to result from 
\0m > being the ground state of H M (in similar fashion to the vanishing of eq. (|5"T|)). 
In other words the total Minkowski energy does not fluctuate and is always equal to 
its eigenvalue zero: vacuum fluctuations carry no energy. Notice also how the ie en- 
codes the above mentioned slight asymmetry between the left and right quadrants: 
(T VV (U < U a ,V)) e is real and positive for V < whereas it is complex and oscillates 
for V > 0. 

In view of the vanishing of the total Minkowski energy (eq. (|141|) ) and of the posi- 
tivity in the region V < 0, the energy in the region V > must integrate to an exactly 
compensating real and negative value. This is not in contradiction with the positivity 
of Rindler energy in the right quadrant, eq. ( |137| ), since the expressions for the Rindler 
and the Minkowski energy differ by the jacobian dv/dV = 1/aV. The oscillations of T vv 
for V > that occur in eq. ( |133| ) as v — ► — oo (which are negligible in the Rindler de- 



scription) are dramatically enhanced by the jacobian in such a way that the Minkowski 
energy in the right quadrant becomes negative, c.f. eq. (|68|) . 

For U > U a , all V, after the atom has made a transition, the Minkowski energy takes 
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the form 



1 g 2 m 2 f jXn x 2 



1 m(iV m + 1) i 7T 
|1 —hx{-aV-ie) 



a 2 V 2 2y / 7rTCo maT 2 maT 2 

\^ e -[\n(-aV -it)} 2 / a?T 2 1 2 1 g - im ln(-aV-ie)/a 1 2 (142) 

It is manifestly real and positive. This is as it should be since we are calculating the mean 
value of the energy in a state that contains one Minkowski quantum. In particular the 
integral fdV (T VV (U > U a ,V)) e is strictly positive (see eq. fl5PP). The time evolution 
has transformed the conditional value eq. ( |140| ) which was complex and carried no 
energy into a real conditional value carrying positive energy. (The change in time of the 
conditional values is further discussed in Part 5). 

Notice how the it prescription in eq (|142|) encodes the asymmetry between the left 
quadrant (proportional to N m + 1) and the right quadrant (proportional to N m ). By 
absorbing the positive Rindler energy m, the two level atom has reduced the negative 
Minkowski energy on the right thereby converting a vacuum fluctuation into a quantum. 
This is summarized in figure 5. 



3.4 The Energy Correlated to a Rindler State 

Up to now we have considered the final state of the atom to isolate certain field con- 
figurations (i.e. those correlated to transitions of the atom). These field configurations 
can also be isolated without making appeal to the accelerated atom through the intro- 
duction of projection operators acting directly on field states, decomposing thereby the 
mean value. To make contact with the previous section we work in Minkowski vacuum 
and consider projections onto states containing certain Rindler quanta. 
We first write unity as 

J = £n, (143) 

i 

where IT are a complete set of projectors. The mean value of Tyy is then decomposed 

as 

(0m\T V v\0m) = (OMl^njVvlOM) 

% 

= ^^ (OAf|njy y |o M ) 



where = (OmIITIOa/) is the probability to be in the eigenspace of IT. As in Section 
Ol the matrix element 



(0 M 


U-iTw 


0m) 


(0 


M 


IT|0m> 



{Tvvhi = /n iTJ tz r (145) 
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is the conditional value of Tyy if the final state is in the eigenspace of Ilj. 

In this section we shall study the properties of some typical conditional values in 
preparation for black hole physics. The physical relevance of these matrix elements will 
be displayed in Part 4 and Part 5. 

We first consider the projector 

H\,R;\,L = a \,R a \,A^RL >< ®RL\a\,Ra\,L (146) 

which projects onto the state containing one pair of rindlerons of Rindler energy A. By 
availing oneself of the identity 

4,r4>,l\Qrl > = —^—4,M a -y,M\°RL > + ~ 6 ( X ~ X ')\°RL > (147) 
a x a x > a x 

it is straightforward to obtain 



<o M | 


iW,l ct>{v)cj>{v')\o M > 


< o M | 


~R\,R;\,l\ 


Om > 



< ORL\ax,Rax,L HV)(I>(V')\0m > 

< 0_RL|ctA,,RaA,I,|OA/ > 

-¥-\m\ V )V\m\ V ) + — — " ( 14 8) 



o^xPx ~~ ' <0 bl \0m> 

It decomposes into two terms. The first depends on the quantum number A and is the 
contribution of the pair of rindlerons selected by the projector ^-x,r-xl- It carries an 
energy density equal to 

lim d v d v ,^-^ XM (VM M (V) = / 2n/ _V 2 (149) 
v^v axpx 27ror (y + ie) A 

The second term is independent of A and appears because except for the mode A, 
Rindler vacuum has been selected (indeed if the projector H-o RL — \®rl >< ®rl\ is used 
only the second term appears). It is convenient to rewrite this term as the sum of the 
expectation value of Tyy in Minkowski vacuum^ plus their difference 



(Tvv) 



RL 



< Orl \ 


Tyy\ 


Om > 


<0rl\ 


Om > 



<0 M \T VV \0 M > + Vjm dydy, / dX -2^ip*_ x>M {V)ipl M {V) 
v-*v Jo ax 



fix 
a x 

< o„|r,„|o M >- T2 {^) ^vTW (150) 



7T / a \ 2 



6 In this section and the following one we shall explicitly write the vacuum expectation value of the 
energy momentum tensor < 0m|?Vv|0m > even though it vanishes. It is kept only to facilitate the 
transcription of these results to the black hole problem where the vacuum expectation of the energy is 
non trivial and must be renormalized carefully. 
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The Rindler interpretation is obtained by considering the Rindler energy density T vv . 
Then eq. (|149|) gives the energy density of the selected rindleron A/27T, the jacobian being 
(dV/dv) 2 = a 2 V 2 . The second term, eq. ( 150 ), is the Rindler vacuum energy p{| which 



is minus the thermal energy density at a temperature a/2ir (Minkowski vacuum contains 
a thermal distribution of rindlerons). The energy in the left quadrant is identical to that 
in the right quadrant since there is a complete symmetry between the two. 

The Minkowski interpretation is completely different. Since the hamiltonian Hm is 
diagonal in u and annihilates Minkowski vacuum Hm\0m >= 0, both eq. (|149|) and eq. 



(|150D contain zero Minkowski energy. Indeed, the pole prescription at the horizon V — 



ensures that their integrals over the entire domain of V vanish as in eq. ( |141| ) . 

The projector Hx,r-,x,l cannot be used in the black hole situation since we have no 
access to the region beyond the horizon nor does it mimick the existence of a detector 
confined to the right Rindler quadrant. Hence we introduce a projector which selects 
the presence of a rindleron of energy A in the right quadrant while tracing over the state 
of the field in the left quadrant: 

n A ,/? = h ®ci{ r \0r >< 0n\a\,R (151) 

where II is the identity operator restricted to the left quadrant and \0r > is Rindler 
vacuum in the right quadrant. The corresponding conditional value of Tyy is given by 



<0 M 
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x,rTvv 


o M > 


< 


M 




M > 



It leads back to eq. ( |148| ) because of the EPR correlations between the two quadrants: 



if there is a rindleron on the right then their necessarily also is a rindleron on the left 
(its partner) with the opposite Rindler energy. This partenaria follows from eq. fl2"5] ) 
where the operators a\ R and a\ L appear in product only. In the black hole problem 
the equivalent EPR correlations will mean that to each outgoing Hawking photon their 
corresponds an ingoing partner on the other side of the horizon. 

An even less restrictive projector specifies only partially the state of the field in the 
right quadrant. One chooses that the final state contains one rindleron on the right in 
the mode A while tracing over all other right rindlerons and over all left rindlerons. The 
resulting projector is 

flx,R = h ® II h',R ® \1\,r >< U,n\ (153) 

X'^X 

where 7a,l is the identity operator restricted to the mode A in the left quadrant and 
|1a,r > is the one particle state restricted to the right mode A. The corresponding 
conditional value of Tyy is 

' + < M \T VV \0 M > (154) 



< M 
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x,rTw 


0m > 
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M 




0m > 



2?ra 2 (V + ie) 2 



The first term is the energy of the rindleron A already obtained in eq. ( |148| ) and eq. ( 152 ). 



The second term is simply the Minkowski vacuum expectation value since no further 
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specification is imposed on the final state. This is why the probability < 0m\^-x,r\0m > 
to be in the eigenspace of IT A r is finite. This is to be opposed to the probabilities 
encountered previously (the denominators of eq. ( |148| ) and eq. ( |152| )) which vanish 
because all the Rindler modes have been specified to be in their Rindler ground state. 
(In physically realistic situations only nonvanishing probabilities will occur. This was 
indeed the case for the accelerated two level atom coupled to the field). 

Nevertheless the conditional values eq. (|152[) and eq. (|154|) are related by a unitary 



relation similar to eq. (|144|) : by taking a set of orthogonal projectors like ITa.^ whose 
combined eigenspace is equal to the eigenspace of IIa^ and summing the corresponding 
conditional values multiplied by the relative probabilities that they occur, eq. ( |154| ) 
is recovered. In order to realize this unitary relation one must select the presence of 
two, three, any number of rindlerons. The corresponding conditional values of Tyy are 
easily obtained and the contribution of each individual selected particle is found to be 
independent (if the particles are orthogonal) of the selection performed on the other 
particles. In other words, for a free field the vacuum fluctuations of orthogonal particles 
are independent of each other. 



We finally consider the selection of a wave packet. Instead of the projector eq. ( |151| ). 
we define: 



n 



vq,\q,R 



h <8> at 



vq,\q,R 



R >< R \a. 



v ,X ,R 



(155) 



where a VOt \ ,R = f^° dXf(X)a\ i R is the destruction operator of a wave packet of right 



rindlerons centered around v = v o and A = Aq. The state IT 



vq,\o,R 



M > is 



n 



1>0,Ao,i? 



M > 



+00 



dXf*(X)a 




+00 



dX' 



P\> ff\i\ t 



\0rl> 



(156) 



where the EPR correlated wave packet in the left Rindler quadrant appears explicitly. 
Note the asymmetry of the wave packets: the induced wave packet in the left quadrant 
contains the factor f3 x ,/ay since it originates from the EPR correlations in eq. (|25|). This 
asymmetry plays a fundamental role when analysing the flux emitted by the accelerated 
detector and the black hole. It is responsible of the fact that Im (T VV (U < U a )) e vanishes 
for V < only, see eq. ( |135| ). The conditional value of Tyy associated with this wave 
packet is 
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v ,Xo,R T VV 


M > 


< 


M 




Om > 
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Om > 



(157) 



The role of /(A) in this equation is very similar to that of c\ in Section [O] (where 
c\ was the Fourier transform of the coupling to the atom, see eq. (0)). However since 
/(A) in eq. ( |157|) contains no negative frequencies, eq. (|157|) is singular on the future 
horizon (see discussion after eq. ([54]) ). 
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4 Black Hole Radiation 



4.1 The Kinematics of the Collapse and the Scattered Modes 



We work in the background metric of a spherically symmetric collapsing star of mass 
M. Outside the star the geometry is described by the Schwarzschild metric 



ds 2 

V, u 
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The specific collapse we consider is produced by a spherically symmetric shell of 
pressureless massless matter. Inside the shell space is flat and the metric reads 



ds 2 = dr 2 - dr 2 - r 2 dQ z 

= dUdv-r 2 dtt 2 
v,U = t ± r 



(159) 



where v is the same coordinate in eq.( |158|) and eq. (|159 ) since on Z~ (u = — oo) space 
time is flat on both sides. The collapsing shell, taken to be thin, follows the geodesic 
v = Vs- The connection between the two metrics is obtained by imposing the continuity 
of r along the shell's trajectory 



dU = du(l — 

Then by choosing v$ = 4M one gets 

du 
u{U) 



2M 



r(u,v s y 



du(l 



AM , 
vs-U' 



(160) 



dU , d , , TT . 
-—{AM-U) 

U — AM In (777) 
y AM' 



(161) 



In the static space time outside the star, the Klein-Gordon equation for a mode of 
the form cp ljTn = -j=L=Y lm (6, tptyfar) reads 



d 2 - e 



'1 - 



2M, 



1(1 + 1) 



+ m A + 



2M 



(162) 



Near the horizon r — 2M << 2M, it becomes the wave equation for a massless field 
in 1 + 1 dimensions. By considering only the s-wave sector of a massless field and 
dropping the residual "quantum potential" 2M(r — 2M)/r 4 the conformal invariance 
holds everywhere, inside as well as outside the star. From now on we shall work in this 
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simplified context and only discuss briefly the differences with the more realistic four 
dimensional case. 

The Heisenberg state is chosen to be the initial vacuum i.e. vacuum with respect to 
the modes which have positive v -frequency on Z~ . These modes are reflected at r = 
and read 

<p U)0fl {v, u) = (e-^ - e-^M) (163) 

Hence, for u > AM (or even on both sides of the horizon for — M < U < M) the state 
of the field tends exponentially quickly (in u) to Unruh vacuum pj, i.e. vacuum with 
respect to the modes 

exp(— iujv) and exp(-^e4M) (164) 

The Schwarzschild w-modes Xx( u ) = e~ lXu /(47rrvA) are needed to analyse the parti- 
cle content of the scattered modes ip^ on X + . In terms of U they take the form 

The exact Bogoljubov coefficients between ip u and xx are given by 

1 fu 



a„, x = (<p u ,xx) = ^jT{l + iAMX)[AM{u; - A)]~ i4MA e ±27rMA (166) 

where the ± is to be understood as + if u > A and — if u> < A. The expression for f3 u> \ is 
obtained by taking A into —A. The asymptotic Bogoljubov coefficients (relating Kruskal 
modes to Schwarzschild modes) are identical to the coeficients relating Minkowski modes 
to Rindler modes (see eq. (plf) et seq.) in the limit uj — > +oo which corresponds 



to resonance at late times u — > +oo (see eq. (|107| )). In this limit the black hole 
emits un correlated (see eq. (^6|) ) quanta at the Hawking temperature l/87rM since 

\Q l n 1 2 _ „-8ttA/A 
\Pu,X/ a w,X\ — e 

Having described the kinematics of the collapse we now turn to the description of 
the energy content of the emitted quanta. The new difficulty lies in the renormalization 
of the energy momentum tensor which must be carried out in curved space times. We 
therefore turn to this point. 



4.2 Matrix Elements in Curved Space-Time 

Wald has proposed a set of eminently reasonable conditions that a renormalized energy 
momentum operator should satisfy By an argument similar to Wald's (or simply 
by verifying that it is in accord with his axioms), it is possible to deduce that T Miy ( ren ) (x) 
can be written in the following way 

T^ T en)(x) = T^(x) - ^(S) {x)I (167) 

where T^ v (x) is the bare energy momentum tensor. The subtraction term t M!/ (s)(^) 
is an (infinite) conserved c-number function only of the geometry at x. It can be 
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understood |36|]|37|] as the (infinite) ground state energy of the "local inertial vacuum": 
that state which most resembles Minkowski vacuum at x. Numerous techniques have 
been developed to calculate t^s) and we refer the reader to [[| for a review. 

In a state, say the Heisenberg vacuum |0 >, the expectation value of T^ v takes the 
form 

< 0|T Mrcn) (x)|0 > = < 0|7>(x)|0 > - V(s)Or) (168) 

where both terms on the r.h.s. are infinite but their difference is finite. 

If in addition one specifies the final state (i.e. on X + ) to be I1|0 >, where II is 
a projector (or more generally the self adjoint operator), the renormalized conditional 
value of Tjjj, reads 

Wu - <0 ,n|o> ^ 169 ) 



Inserting eq. (|167|) into this expression yields 



(T^{x))u = < Q |ri|Q > V(S)W (170) 

Then by expressing T^ u (x) in terms of the operators which annihilate the Heisenberg 
vacuum one obtains 

(T»v(x)) u = fdu / o ^ <0 ^^ 0> f^(x) [<p*M + < 0\T^ {len) (x)\0 > (171) 

where T^x) is the classical differential operator which acting on the waves eq. 
(|163|), gives their energy density. 



The renormalized conditional value contains two contributions. The first term, the 
fluctuating part, depends on the particle content of the state specified by II. Contrari- 
wise, the second one is the mean energy density of the Heisenberg vacuum eq. (|168|) 
obtained when no specification on the final state is added. 

The formula eq. ( |170| ) warrants a few additional comments. First notice that their 
are parts of (T /J ^) n that are entirely contained in the subtraction. Most notably there 
is the trace anomaly and those components of the energy momentum tensor which are 
related to it by energy conservation (in two dimensions they are T UUyV and T VVyU ). These 
parts are independent of n or, expressed differently, do not fluctuate. 

An additional (and related) feature concerns the absence of correlations between T uu 
and T vv . Not only shall this give rise to the particular structure of vacuum fluctuations 
that extend back to Z~, but it also implies that on the horizon the in-going flow and 
the out-going flow fluctuate independently (for instance the specification of an outgoing 
particle on X + does not affect T vv outside the star and in particular on the horizon 
r = 2M). This last effect disappears partially when considering the potential barrier 
that occurs in the wave equation eq. 



4.3 The Conditional Value of the Energy Density 

In the absence of specification of the final state, < 0|T^( ren )|0 > describes the mean en- 
ergy content carried by Hawking quanta. We remind the reader that < 0|T^,( ren )(x)|0 > 
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is regular on the future horizon U — and that Hawking radiation can be conceived as 
the matter response that gives regular mean energy densities on the horizon. We refer 
to ref. || for further discussion of the mean flux. 

In order to describe the fluctuations around the mean flux, one inquires into the 
conditional energy density when the final state is a particular out-state which arises in 
the rewriting of the Heisenberg state |0 > into states with definite energy A on Z + , see 
eq. fl2"o]). However in the collapsing geometry, an external observer does not have access 
to the region of space time beyond the horizon, hence the specifications of the final state 
that he can perform are restricted to an incomplete (U < 0) region of space time and 
are therefore incomplete as well. 

A possible specification is, for instance to use the projector, see eq. ( |151| ). 

n A = I L <g> a\\B R >< B R \a x (172) 

where II is the identity operator restricted to the inaccessible region U > and \B R > 
is Boulware vacuum in the region U < 0. a\ creates a Schwarzschild outgoing photon. 
The corresponding conditional energy reads, see eqs. ( |148| ), ( |150| ) and ( |171| ), 



+ 



< B\T U JX)\0 > \m / \ I 

1 -<0|T^(x)|0> 



< BO > 



< 0|T Mrcn) (x)|0 > (173) 



where the modes <fx,K are defined from the modes e lU3U / V 'Ahuj , eq. (|163|), as the 



Minkowski modes eq. ( 15|) are defined from the Minkowski plane waves eq. (|6|). 

The first two terms arise from the specification of the final state whereas the third 
term is the mean energy. The first term is equal to the energy of the photon A. The 
second one is the difference of energy between Boulware vacuum and the Heisenberg 
vacuum. This term appears, as in eq. ( |150| ), because one has specified that, apart from 
A, their is no other photon emitted. This is why this term is singular on the horizon and 
why the probability to be in the eigenstate of Ha vanishes in the absence of backreaction 
(in the semiclassical approximation, it is of order e _A/2 where M 2 is approximately the 
total number of photons emitted). 

A more reasonable specification because it has a finite probability of occurring con- 
sists in tracing over all the photons except the photon A which is imposed to be present 
(in the Rindler problem this corresponds to the projector eq. ( |153| )). Then the condi- 
tional energy is simply, see eqs. ( |154j ) and (|171|) , 

(^(z))n A = ^^{x) [pIkV-x,k] + < 0|T Mren) (x)|0 > (174) 



'A 



Having traced over all the other photons, the second term of eq. (|173|) is absent 



in eq. ( |174p . Nevertheless this term can also be constructed as the sum of conditional 



values that specify completely the state times the probability that they occur (in similar 
manner to the the unitarity relation eq. ( |144| )). In this way the difference of energy 
between the Heisenberg vacuum |0 > and Boulware vacuum \B > is realized as the sum 
over all possible radiated photons times the thermal probabilities that they occur. 
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Finally we consider an inertial two level atom at large distance from the black hole. 
The specification of the state of the radiation is carried out indirectly by requiring, as 
in Section 3.2, that the atom get excited. In this case also the final radiation state is 
partially specified, since the detector is coupled to a finite set of modes. One finds that 
the conditional energy contains again two terms 

(T^) n+ = (T^) e + < 0|T^ (ren) |0 > (175) 

where the fluctuating part, in terms of the 4> m operators (eq. (^)) is, see eq. ( |127] ), 

( T ^)e = °l4An : T, u : |0 > (176) 

e 

where : T^ v : is the energy momentum operator normal ordered with respect to the in 
operators defining Heisenberg vacuum |0 >. We shall display the properties of (T^) in 
the next section. 

We note already that the specification of the final radiation state by the correlations 
to a transition of the two level atom gives rise to finite energy densities on the horizon 
only if the coupling to the field decreases faster than e ~ u ^ M : c.f. discussion after eq. 
(|157|) and (|54]). Therefore, the specification of a mode (eq. ( |174j )) or a wave packet made 



out of positive A frequencies only gives rise inevitably to singular energy densities on the 
horizon. 

It is also interesting to speculate about the nature of the in-going vacuum fluctua- 
tions when they cross the future horizon U = 0. They could be analysed by selecting 
the presence of ingoing quanta near the horizon. A "natural" set of modes to select near 
the horizon are Kruskal f-modes. One is therefore led to consider the Kruskal vacuum 
fluctuations in Schwarzschild vacuum, which is similar to considering Minkowski fluctua- 
tions in Rindler vacuum. If space time were the full Schwarzschild manifold, these would 
present a singularity on the past horizon that could be smoothed out using wave packets. 
Since space time is not the full Schwarzschild manifold (there is no past horizon) the 
star's surface will play the role of past horizon and one expects large energy densities in 
the outermost layers of the star. 

4.4 From Vacuum Fluctuations to Black Hole Radiation 

We now turn to the the fluctuating part of the conditional energy correlated to the 



transition of the two level atom, i.e. the term {T^ v ) e of eq. ( |175|) . 

For s-waves, when one neglects the residual potential of the d'Alembertian eq. 
(T /ii/ ) e is completely independent of the geometry. Furthermore, the conformal invariance 
of the field makes the mapping of the results obtained in the Rindler problem to the 
present problem straightforward. Let us choose therefor the time dependent coupling 
fit) of our detector given by the gaussian switch off of Section 3.2, see eq. (|57D and eq. 
(|63|). More precisely, the coupling is such that the atom, of resonance frequency A = m 
will be excited around the retarded time u = uq. Hence its Fourier components are, see 
eq. © 

Cx = Jdt f(t)e- vmt e iXt = D—e iXuo e-^- m)2T2 / 2 (l - e - 2nX/a ) (177) 
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The spread in time is At = Au = T and Uq is taken well inside the region u > 0, see eq. 
(|161|) , where the isomorphism of the scattered waves and the Kruskal modes is achieved. 



If the two level atom is found excited after the switch off, due to the vanishing of 
the modes at r = 0, see eq. (|163|) , and the light-like character of the propagation, the 



correlated radiation state eschews from a spherically symmetric vacuum fluctuation on 
X~ . This fluctuation is located in a region 



V - V r , 



\AU\ = \Aue- uo/4M \ ~ Te~ uo/4M (178) 



where v = v ^ (= in our collapse) is the light ray that shall become the future horizon 
U — 0. Indeed this localization is furnished by the v dependence of the conditional 
energy density on X~ which reads (see eqs. ( |176| ) and (|140| )) 



g 2 m 2 



{T vv (l ,v)) e = —— d\ d\'c* x c x >- — 7 =Jh X '(n x + 




1 16M 2 m . 
W^V?T (iVm+1) 6XP 



AM n ,-v-ie. , n2 
(ln( AiLf )+uq) 



T AM 



(179) 



Where the width of the gaussian factor gives eq. ( |178j ). As in eq. ( |140| ), the e specification 



of the log ensures that the total energy carried by this fluctuation vanishes, see eq. 
( 14ip . Thus we see that the analysis of the fluctuations by an inertial observer near 



X - is isomorphic with what was called the Minkowski interpretation in Part 4. As in 
the accelerated case, the energy density is enhanced by the jacobian du/dU = e u ^ iM 
centered around u — Uq which appears here as 1/v 2 when the reflection at r = is 
taken into account. Hence after a w-time of the order of 4MlnM, the energy density 
in T vv (rescaled by 4irr 2 ) become "transplanckian" and located within a distance Av 
much smaller than the Planck length (If one does not rescale the transplanckian 
energies only exist in a region of finite r which nevertheless increases exponentially with 
«o). The dramatic consequences that these transplanckian energies might introduce 



are discussed in ref. ||19|| . In that article it is argued that the nonlinearity of general 
relativity cannot accommodate these densities and that a taming mechanism must exist 
if Hawking radiation does exist. 

After issuing from X~, the vacuum fluctuation contracts until it reaches r = and 
then reexpands along U = const lines. Upon crossing the surface of the star in a region 
AU ~ y e -«o/4M cen tered on the horizon, it separates into a piece (the partner) that falls 
into the singularity, carrying a negative Schwarzschild energy equal to — m(l + N m ) (see 
eqs. |135| and |139j ), and a piece carrying positive energy equal to m(l + N m ) that keeps 
expanding and escapes to X + to constitute the quantum that induces the transition of 



the atom (see eq. |133| and figure 4). The analysis performed by an inertial asymptotic 



observer near the detector, on X , is isomorphic with the Rindler interpretation of 
Section 3.3 since the gravitational red shift replace exactly the role of the Doppler 
accelerated one in the accelerated place. For instance, one finds readily that the total 
energy carried by the fluctuating (T uu (l + , u)) e is indeed m(l + N m ) as in eq. ( |137| ). 

Similarly, if the two level atom is found in its ground state after the switch off, its 
wave function is correlated to the absence of the Hawking photon specified by c x . In 
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that case, one would find near X~ a vacuum fluctuation whose energy content is exactly 
the opposite of the previously considered case (times P g /P e ). Near X + it would therefore 
contain a negative energy flux of total energy —m(N m + l)P e / (1 — P e ) encoding the fact 
that their are quanta absent from the thermal flux emitted by the black hole. 

If more realistically, we take a two-level atom coupled locally to the field (i.e. coupled 
to all the modes I > 0), it will select particles coming out of the black hole in its direction. 
Then the picture that emerges is essentially the same as for an s-wave except that on Z~ 
the vacuum fluctuation is localized on the antipodal point of the detector. The created 
quantum and its partner, are on the same side and not antipodal (with respect to each 
other) because they have opposite energy. 

We now turn to the description in the intermediate regions in order to interpolate 
between the descriptions between I~ and I + . One possible interpolation consists in 
using a set of static observers at constant r. Then the "Rindler" description would be 
used everywhere outside the star. However a difficulty arises in this scheme if one really 



considers a set of material "fiducial" [20] detectors at constant r. For upon interacting 



with the field and thermalizing at the local temperature ^J 2 m) 8^m ^ ne detectors will 
emit large amounts of ultraviolet Kruskal "real" quanta (see Section 3.3 wherein it is 
shown how the accelerated atom transforms vacuum fluctuations into "real" quanta). 
The backreaction of these on mass shell quanta cannot be neglected and, as already 
stated, cannot be evaluated owing to the transplanckian energy they carry. 

An alternative interpolation consists in giving the value of T^ v in the local inertial 
coordinate system (Riemann normal coordinates). This stems from the idea that local 
physics should be describe locally in such a coordinate system. This approach has 
been used in defining the subtraction necessary to renormalize the energy momentum 



tensor ||36|| | S7| . In the two dimensional model the local inertial coordinates are easy 
to construct. Since u = r(u, v) is an affine parameter along the geodesies v = constant, 
a natural way to represent the outgoing flux outside the star is as 

Tuu(u) = ( ^ U ^ V) ") Tuu(u(r,v)) (180) 

This is represented in a Penrose diagram in figure 6 and Eddington-Finkelstein coor- 
dinates in figure 7. The inertial coordinate u will come up very naturally in the next 
section. 



4.5 The Gravitational Back Reaction 

Up to now in this Part we have presented the properties of the conditional values of 
Tfja,. We now investigate how these matrix elements intervene in physical processes. An 
example of the role of these matrix elements was given in section |3.2| . Here we shall 
specifically discuss gravitational back reaction effects to black hole radiation. 

A simple way to understand the role of the conditional values of T^ v is to imagine a 
change in the background geometry g^ v — > g^ u + 5g^ u . This change modifies both the 
mean values of the flux as well as more detailed properties such as the probability to 
find a specific photon on 1 + . We focus on this later change. 
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The new probability can be computed, in the interacting picture, around the back- 
ground with the Hamiltonian given by 

J d 4 xH mt = ~J d 4 x^5g^T^ (181) 
The new probability is then, see eq. ( |144[ ), 

Pg + s g = (Ole^^^ne-^^'^IO) (182) 

where II specifies the state of the radiation field on 1 + . To first order in 5g tiV the relative 
change in probability is 

p 9 +s g -Pg _ (o\n(-ifd 4 xH int )\ty 



p g (o|n|o) 

= J d 4 xy^5g^ lm[(TM (183) 

It is thus the imaginary part of ( T Miy )n only which controls the change in probability 
induced by 5g iiV . Furthermore, since the background part of conditional energies (the 
second term of eq. ( 171 )) is by construction real, only the fluctuating part, which 
depends explicitly of the selected quantum, contributes to P g +s g — P g - 

To illustrate how the various properties of the fluctuating part of (T^n intervene in 
such an expression, let us take the simple example wherein 5g iiv is due to the infall of an 
additional light like shell of mass 5m at time v = v ' with v' > v$- Then for v$ < v < v' 
the metric, eq. (|158|) and eq. ( |161| )), is unchanged 

d s 2 = (l - ^) dv 2 - Idvdr - r 2 d 2 n (184) 

whereas for v > v' it is 

ds 2 = il \dv 2 -2dvdr~r 2 d 2 Q (185) 

where we have used for obvious convenience the Eddington Finkelstein coordinates v 
and r. 

The change in the action S = J d 4 x y /—g^g' JlU d fl (f)d u (j) is for s- waves 

//•+oo /*oo 
d 4 xH int = dv dr Aitr 2 {5m I r)d r (j)d r (j) 

Jv' JO 

/*+oo /*oo 

= 5m / dv dr AnrTuu(r,v) (186) 
Jv' Jo 

where is given by eq. (|180|). As emphasized at the end of section PO it is the 



energy momentum in Riemann normal coordinates which appears automatically in such 
problems since the response to a local change in the geometry is local as well. 

In order to compute the change in probability due to Hint, we first recall that the 
imaginary part of (Tqz)ti vanishes on the other side of the horizon, for r < 2M (see 
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discussion after eq. ( |139| )). From eq. ( |186| ) we understand that this is dictated by 
causality: a change in g^ v in the region r < 2M cannot affect the probability to find a 
specific Hawking photon on this side. 

In the case when v' = vs, one simply has a shell of mass M + 5m. In this case 
we know exactly the change in probability since the probability of finding a photon of 
frequency A is 

PM + s m = |/WI7K,a| 4 = e~ 8 ^ M+5m \l - e - 8 ^ M+s ^) (187) 

hence 

5P (l-2e~ 8nXM ) 

T = - 8nX6m (l - e^AM) ( 188 ) 

Thus in this case the transplanckian character of the energy density is washed out by 
the integration in eq. (|186|) . This can be verified explicitly by evaluating the integral 
near r = 2M and v = v'. There, one can replace dr by —dU/2 whereupon by making 
appeal to the vanishing of the integrals JdU(Tuu)~(i and Jdulm(T uu )ji one finds that the 
integral at fixed v J<ir47rrIm(T fifi )n does not scale like e U//4M . 

Furthermore at large r the integral, eq. (|183|) , vanishes once more since J du\m{T uu )n 
vanishes. Hence if the additional mass crosses the photon trajectory when it is on mass 
shell at r >> AM there is no modification of the probability of creating the photon. 

We have thus obtained a local description of the quantum matter response to a mod- 
ification of the classical background geometry. To address the quantum gravitational 
back reaction to the creation of a Hawking photon one should treat Sg^ as a quantum 
operator [IT] [II]. For spherically symmetric radiation and spherical symmetric gravita- 



tional fields the relation between Sg^ and T^ v is a constraint, i.e. Sg^ is completely fixed 
by T M ;,. One can then envisage the backreaction as an iterative scheme which ultimately 
should be treated self consistently. 

The first step in this procedure is very simple. It consists in taking bg^ v to be the 
"position coordinate" of the additional system (the oscillator) introduced in Section 



2.5. Then, as for the oscillator, see eq. (|132|) , the mean value of 5g, w is obtained by 



integrating Einstein's equations with the mean energy momentum tensor as a source. 



This corresponds to the linear approximation to the semiclassical solution |fLq| . But one 
can also evaluate the "mean" conditional value of 5g^ u . This conditional change in the 
metric is obtained by integrating Einstein's equations with (T fJiU ) Tl as source. [] Since the 
total energy carried by the conditional value of vanishes from X~ till the emergence 
of the fluctuation from the star after reflection on r = 0, {5g^ v ) n) the conditional value 
of Sg^ u , will vanish outside the interval Av eq. |178| centered around vs- Within that 
interval the precise shape of {5g^ u ) n will depend on the particular choice of selected wave 
packet by the projector n. On the contrary, outside the star, for r > 4M and u > u 
(i.e. in the middle of the two members of the pair), (5g^ u ) u will encode the mass loss u 
and in fact describes a new classical (real valued) Schwarzschild space where the mass 
is M - to. 



7 See however the different ways in which Re(T uu )n and Im(T utt )n enter in eq. ( |f 32| ). This might 
lead to interesting effects. 
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The next step consists in taking into account the effect of Sg^ on the production 
of the Hawking photon itself. This gravitational self interaction can be encoded in 
a interaction hamiltonian of the form Hi nt = T^D^^Tap where D is the linearized 
gravitation propagator. To calculate 5P/ P due to this self interaction one must confront 
the infinities which arise in matter loops. This will not be done here. 

However there are some simple question wherein H int does come in which do not 
involve loops. One such question is the effect of the creation of a first Hawking photon 
on subsequent ones. Suppose one calculates Pn Xl u . x „ 2 where ^-\ 1} urM,u2 specifies the 
presence on X + of a photon of frequency Ai located near u\ and another photon of 
frequency A2 located near u-i- Then to first order in Hi nt , the relative change in the 
probability to find the two photons is 

$ p n x x 

Pu i -Pn \ 

where, in the absence of gravitational coupling, the probability of finding two photons 
factorizes into the individual probabilities. In this expression once more there are infinite 
loops. However upon taking derivatives the quantity (d/d\i)(d/d\2)(SP/P) is finite. 
The effects of particle Ai, u± on particle A2, «2 are isolated from other effects. This finite 
quantity can then be expressed in terms of the products of the one-particle conditional 
values 

(^(^))n Al , ul ^ a/3 (^,2/)(T a/3 (y))n A2 , tt2 

(^(^^(y))^^^^!/)^^)^^)^ (190) 

In conclusion we have shown that the conditional values of T uu enter into tree graphs. 
In order to understand the role of the transplanckian frequencies in Hawking radiation 
and how they are tamed by quantum gravity one should confront loops and the infinities 
they involve. We hope to report on this in subsequent work. 



< o|n 


Ai ,«i ; A2 ,1x2 Hint 


> 


^ 0|n.\i,Ui;A2,M2 


> 



(189) 



5 Post Selection, Weak Measurement 
5.1 Introduction 

In section |3.2| we showed succinctly how the conditional values of Tyy control the first 
order perturbation onto an additional system. The aim of this section is to present a 
self contained discussion devoted this result. We shall work in complete generality and 
make no explicit reference to the accelerated system nor to black hole radiation. 

This analysis of the conditional values (i.e. non diagonal matrix element of an oper- 
ator) was first carried out by Aharonov et al. |IIJ in the context of measurement theory. 
In essence they studied the first order backreaction onto an additional system which 
they took to be a measuring device. But the formalism is more general. In the case of 
pair production the additional system could be the external electric or gravitational field 
which is now described quantum mechanically. Moreover this formalism can be used to 
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study the self interaction of the pairs without introducing the additional system. This is 
because, when the first order (or weak) approximation is valid, the backreaction takes a 
simple and universal form governed by a c-number, the conditional value of the operator 
which controls the interaction (called by Aharonov et al. the "weak value"). 

In section 5.2 we implement the specification of the final state(s) in a rather formal 
way by acting with projection operators which select the desired final state(s). Aharonov 
et al. call this specification of the final state, a "post- select ion" . 

In section 5.3 we show how post-selection may be realized operationally following 
the rules of quantum mechanics by coupling the system to be studied (S) an additional 
system in a metastable state (the "post-selector" PS) which will make a transition only 
if the system is in the required final state(s). The conditional value of an operator 
obtained in this manner changes as time goes by from an off diagonal matrix element to 
an expectation value, thereby making contact with more familiar physics. This extended 
formalism finds important application when considering the physics of the accelerated 
detector since the accelerated detector itself plays the role of post selector. 



5.2 Conditional (or Weak) Values 



The approach developed by Aharonov et al.||r2| for studying pre- and post-selected en- 
sembles consists in performing at an intermediate time a "weak measurement" on S. 

The system to be studied (S) is in the state \ipi > at time ti (or more generally is 
described by a density matrix p,;). The unperturbed time evolution of this pre-selected 
state can be described by the following density matrix 

Ps(t) = Us(jt,U)\if>i >< if>i\U s (ti,t) (191) 

where Us = exp(-iHst) is the time evolution operator for the system S. The post- 
selection at time tf consists in specifying that the system belongs to a certain subspace, 
7i s , of Tis- Then the probability to find the system in this subspace at time tf is 

P n o = Tr s [n s p{t f )} = Tr s [n s U s (t f ,t t M >< ^\U s (t t ,t f )] (192) 

where 11^ is the projection operator onto H s and Tr s is the trace over the states of 
system 5*. In the special cases wherein the specification of the final state is to be in 
a pure state \ipf > (i.e. H| = \ipf >< ipf\) then the probability is simply given by the 
overlap 

P f = \< ^f\Ug(t f ,UM >l 2 (193) 
Following Aharonov et al. we introduce an additional system, called the "weak 
detector" (WD), coupled to S. The interaction hamiltonian between S and WD is 
taken to be of the form H s _ WD (t) = ef(t)A s B WD where e is a coupling constant, f(t) 
is a c-number function, As and B>wd are hermitian operators acting on S and WD 
respectively. 

Then to first order in e (the coupling is weak), the evolution of the coupled system 
S and WD is given by 

p(t f ) = \V(t f )><*(t f )\ 
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where 



*(*/) > = [Usitf^UwDitf^-it f'dt U s (t f ,t)U WD (t f ,t)f(t)A s B WD x 



U s (t,U)U WD (t,U) \A >\WD > 



(194) 



where Us and Uwd are the free evolution operators for S and WD and \WD > is the 
initial state of WD. Upon post-selecting at t — tf that S belongs to the subspace Ti, s 
and tracing over the states of the system S, the reduced density matrix describing the 
WD is obtained 



PwD(t f ) = Tr s 
To first order in e, it takes a very simple form 



PWD 



(195) 



where 

\VwD(tf) > 



U W D(tf,ti) -ie dt U WD (tf,t)f(t){A s (t)) Il oBwDUwD(t,t i ) 



\WD > 
(196) 



where P n o is the probability to be in subspace 7i s and 



Tr s 


n s Us(tf,t)A s Us(t,t i )\^ i >< ^\u s (t 




Tr s 


'iPsUsfauMi ><^i\Us(t u t f )_ 





(197) 



is a c-number called the weak or conditional value of A. If one specifies completely the 
final state, Ii s = \ipf >< ipf\ then the result of Aharonov et al. obtains: 



<^s(*)> 



<jj f \Us(t f ,t)A s U s (t,t i )\ij i > 
< ipflUsitf,^)^ > 



(198) 



The principal feature of the above formalism is its independence on the internal structure 
of the WD. The first order backreaction of S onto WD is universal: it is always 
controlled by the c-number {A s (t)) u , the "weak value of A\ Therefore if S is coupled 
to itself by an interaction hamiltonian, the backreaction will be controlled by the weak 
value of H int in first order perturbation theory. For instance the modification of the 
probability that the final state belongs to H s is given by the imaginary part of (iJ int ) n o . 
Indeed 



U° s (1 -ijdt H int )p t {l + i Jdt Hi^) 
2 Im(tf int ) n o) 



(199) 



The weak value of A is complex. By performing a series of measurements on WD and 
by varying the coupling function f(t), the real and imaginary part of (As(t)) n could in 
principle be determined. Here the word " measurement" must be understood in its usual 
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quantum sense: the average over repeated realizations of the same situation. This means 
that the weak value of As should also be understood as an average. The fluctuations 
around (As(t)) n are encoded in the second order terms of eq. ( |194| ) which have been 
neglected in eq. ( |194j ). 

To illustrate the role of the real and imaginary parts of (As) u , we recall the example 
of Aharonov et al. consisting of a weak detector which has one degree of freedom 
q, with a gaussian initial state < q\WD > = e~ q / 2A ,— oo < q < +00 (see also the 
example of Section 3.2). The unperturbed hamiltonian of WD is taken to vanish (hence 
UwD(ti,t 2 ) = 1) and the interaction hamiltonian is Hs-wDif) = e$(t — t)pA$ where p 
is the momentum conjugate to q. Then after the post-selection the density matrix of 
the WD is given by Pn\WD(tf) >< WD(tf)\, see eq. ( |196| ), where, to first order in e, 
\WD(t f ) > is 

\WD(t f )> = (l-iep(A s (t)) u )\WD> (200) 
whereupon the conditional value of q and p are 

(q) u = eRe(A s (t)) n 

(p) n = elm(A s (t)) u /A 2 (201) 

Thus the real part of (As(t)) u induces a translation of the center of the gaussian, the 
imaginary part a change in the momentum. Their effect on the WD is therefore mea- 
surable. The validity of the first order approximation requires e\(A s (t)) u \/ A << 1. 

It is instructive to see how unitarity is realized in the above formalism. Take W s to 
be a complete orthogonal set of projectors acting on the Hilbert space of S. Denote by 
Pj the probability that the final state of the system belong to the subspace spanned by 
1T5. and by (As(t)) u j the corresponding weak value of A. Then the mean value of As is 

< ^\A s (tM > =J2 P Ms(t))ni ( 202 ) 
3 

Thus the mean backreaction if no post- select ion is performed is the average over the post- 
selected backreactions (in the linear response approximation). Notice that the imaginary 
parts of the weak values necessarily cancel since the l.h.s. of eq. ( |202| ) is real. Equation 
(|202|) is the short cut used in the main text to obtain with minimum effort the weak 
values. 

5.3 Physical Implementation of Post Selection 

Up to now the postselection has been implemented by inserting by hand the projector Tl s . 
Such a projection may be realized operationally by introducing an additional quantum 
system, a "post-selector" (PS), coupled in such a way that it will make a transition 
if and only if the system 5* is in the required final state. Then by considering only 
that subspace of the final states in which PS has made the transition, a post-selected 
state is specified. This quantum description of the post-selection is similar in spirit to 
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the measurement theory developed in ref. |]13fl : by introducing explicitly the measuring 
device in the hamiltonian the collapse of the wave function ceases to be a necessary 
concomitant of measurement theory. 

We shall consider the very simple model of a PS having two states, initially in the 
ground state, and coupled to the system by an interaction of the form 



H S -ps = \g{t){a*Qs + aQ\ 



(203) 



where A is a coupling constant, g(t) a time dependent function, at the operator that 
induce transitions from the ground state to the exited state of the PS, Qs an operator 
acting on the system S. The postselection is performed at t = tf and consists in finding 
the PS in the exited state. 

For simplicity we shall work to second order in A (although in principle the interaction 
of PS with S need not be weak). In the interaction representation, the wave function 
of the combined system S + WD + PS is to order e and to order A 2 

Te -ijdt H s ^ WD {t)+H S -P S (t)^. > \ WD >| 0p5 > = 

l-ijdt (H S - W D(t) + H s -Ps(t)) 

~\l dt l dt ' T l H s-Ps(t)H S -Ps(t')} - J dtj dt'T[Hs-wD(t)Hs-Ps(t')} 

+ \\ dt \ dt ' J dt"T[Hs-Ps(t)Hs-Ps(t')Hs-wD(t")]\ >\WD >\0 PS >(204) 

where \0ps > is the ground state of PS and T is the time ordering operator. The 
probability of finding the PS in the excited state at t — tf is, at order A 2 , 

Pe = A 2 < A\ J dtg(t)Ql J dt'g{t')Q s \A > (205) 

Upon imposing that the PS be in its excited state att — tf the resulting wave function 
is, to order e and A 2 , 



-i J dt\g(t)Q s (t) 
dt f dt'T[ef(t)A s (t)B WD (t)\g(t')Q s (t')} 



\if>i >\WD >a t |0 P5 > (206) 



Making a density matrix out of the state ( 206 ), tracing over the states of S and PS 



yields the reduced density matrix pwD = P e \^WD >< ^wd\ of WD, where, to order e 
\$>wd > is 





r r tf 


\^WD > = 


l-ie 




Jti 



dif(t)B WD {t)(A s (t)) ( 



\WD > 



and 



(As(t)) ( 



<i> i \ldtg(t)QUt)fdt'g(t')T A s (t)Q s (t>) \^> 



<^\J dtg{t)Q\{t) J dt>g{t>)Q s {t>) |^ > 



(207) 



(208) 
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Note how the weak value of As results from the quantum mechanical interference of the 
two terms in eq. (|206|) . 

There are several important cases when the time ordering in eq. ( |208| ) simplifies. If 
g(t) is non vanishing only after t = t, ie. WD interacts with S before S interacts with 
PS, then (As(t)) e takes a typical (for a weak value) asymmetric form 

(A = <HI dtg(t)Q s (t) J dt'g(t')Q s (t')A s (tM > 
<^\Jdtg(t)Ql(t)Jdt'g(t')Q s (t')\^> 

If in addition g{t) = S(t — tf) and Qs = II^, eq. ( |197| ) is recovered using (IT5.) 2 = Yl s and 
tf > i. This is expected since in this case the post-selector has simply gotten correlated 
to the system in the subspace 7i s 

If on the other hand g(t) is non vanishing only before t = t, ie. S interacts with PS 
before WD interacts with S, then the time ordering operator becomes trivial once more 
and eq. ( [208|) takes the form 

(A = < HIdtg(t)Qlms(t)Idt'g(t')Q s (t'M > 

< iPi\Jdtg(t)QUt)fdt>g(t>)Qs(t>M > 

This is by construction the expectation value of As if the PS has made a transition. It 
is necessarily real contrary to eq. ( [209| ). 

Finally, the weak value of As if the PS has not made a transition can also be 
computed. Once more the two cases discussed in eqs ( [209| ) and ( |210| ) are particularly 
simple: if g(t) is non vanishing only after t = t one finds 

(Mt)) d = T ^y(<HMA> 



-A 2 Re 



< Vi| / dtg(t)Q s (t) J dag{t')Q s {t')As{t)\^ 



> 



(211) 



On the other hand if g(t) is non vanishing only before t = t one finds 

(Mt)) d = <^l4(i)li> 



--A 2 Re 
2 



< ^\A s (t) ( dt ( dt'Tg(t)Ql(t)g(t')Q s (t'M > 



(212) 



These are related to the mean value of As and to eq. (p08|) through the unitary relation 
eq. ( [202j ): if g(t) is non vanishing only after t = t 

P e (A s (t)) e + (1 - P e )(A s (t)) d = < ^i\A s (t)\^ > (213) 
and if g(t) is non vanishing only before t = t 

P e {A s (t)) e + (l-P e ){A s (t)) d = <^\e ljdt Hs - ps A s {t)e- i I dt Hs ~ PS \^> = 



< A\A s (t)\A > + A 2 < A\ J dtg(t)Q s (t)A s (i) J dt' 'g(t)QWM > 

-X 2 Re < ipi\A s {t) JdtJ dt'Tg(t)QUt)g(t')QUt'M > (214) 
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where the r.h.s is the average value of As before (eq. ( |213| )) and after (eq. ( |214| )) the 
detector has interacted with S. 

Acknowledgements. The authors would like to thank R. Brout, F. Englert, S. Popescu 
and Ph. Spindel for very helpful discussions. 



6 Appendix: The V Term 

We recall the this term arises from the following decomposition of the second g 2 Born 
term in \ip^(t = +00) > = e^I dtdxHint \0 M >\- > 

r+oo r-T 

-g 2 m 2 / dr / rfr7(r) e - imT 0(r)/*(r / )e + ^ 0(r')|O M >|- > 



g 2 m 2 f+°° 



2 



00 J —oo 
oo 



00 



dr / dr'f(T)r(r')e- m ^ 'V(r)0(r') [1 + e(r - r')} \0 M >\- > 



2 

where 



Gf 77*1 

^r-4L4> m \OM >|- > - g 2 m 2 V\0 M >\-> (215) 



r+oo r+oo 

v = - drj dnf{T 2 )r{n)e{T 2 - n)e-^-^<f>{T 2 )4>{n) (216) 

Z J —oo J —oo 

and where e(r 2 — T\) = 9(r 2 — t x ) — 6{ji — r 2 ). 

To explicitize the role of the V term, it is appropriate to compute the energy density 
carried by it when the initial state is |0jy >|— >. One finds, to order g 2 , 

(T VV (V)) V = -g 2 m 2 Re[<0 M \T vv (V)V\0 M >] 
= -g 2 m 2 < M \ [T VV (V),V]_ \0 M > 
n^m^ f+°° r+oo 

= - 9 — / drj dnfMrWefo-Tje™**^ 

/ J —oo J —oo 

< A /| [T vv (V),<P(r 2 )<f ) (r 1 )]_ \0 M > (217) 

where we have used the antihermitian property of T>: T>^ = —T>. (Tyy{V)) v enjoys the 
following properties. 

1. Being a commutator, {T vv {y)) v is causal (see eq. ([49])), and vanishes in the left 
quadrant V < contrary to (Tyv{V)) e and (Tyv(V)) . 

2. (Tyy{V)) v carries no Minkowski energy since the hamiltonian Hm, eq. (f|), 
annihilates Minkowski vacuum. 

3. {T V y{y)) v carries no Rindler energy since H R (the boost generator given in eq. 



([23])) annihilates Minkowski vacuum. Therefore by virtue of 1., the Rindler energy in 



the right quadrant (V > 0) vanishes 

r+oo 

dv (T vv {v)) v = (218) 
Thus (T vv {y)) v is, at most, an energy density repartition. 
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4. When /(r) = 1 for all r, (T vv {v)) v vanishes identically. To prove this one evaluates 
the commutator in eq. ( |217|) and one finds 



{T w {v)) v = 2g 2 m 2 / dr 2 e(r 2 - v)Re f {r 2 ) f* {v)e~ vm{r ^ v) < O M |0(r 2 )z9^(i;)|O M > 

J -co J 

(219) 

where we have used the commutation relation 

[T vv {v), 0(r 2 )0(r 1 )]_ = -2i5{n - v)(f>{r 2 )d v (f>{v) - 2i6(t 2 - v)d v 4>{v)(j>{T l ) (220) 

and the antisymmetric character of e{r 2 — T\). Since the expectation value < Om|0(t"2)<K^)|O.a,/ > 
is evaluated along the accelerated trajectory eq. (p6|) , it is a function of r 2 —v only. There- 
fore the integrand of eq. ( j219| ) is an odd function of r 2 — v and the integral vanishes. 
Hence {T vv (v)) v is an energy repartition which is concerned only with the transients 
induced by the switch on and off effects. 

5. When /(r) is a slowly varying function with respect to both 1/m and 1/a (c.f. 
the discussion associated with eq. (|59|), (T vv (y)) v is smaller than the contribution of 
Re[(T OT (v)0 m 0j n )] by a factor 1/aT except near the edges of the interaction period where 
/(t) almost vanishes. This can be seen by developing f(r 2 ) given in eq. |219j in a series 
around t 2 = v and evaluating the magnitude of the first non vanishing term, i.e. one 
treats the variations of the switch off function /(r) as an adiabatic effect. One finds 
that indeed the V is smaller than Re[(T tw ( , u)</> m </>] n )] except when r > aT 2 . 

References 

[1] SW. Hawking, Nature 248, 30 (1974) 
Commun. math. Phys. 43, 199 (1975). 

[2] W.G. Unruh, Phys. Rev. D14 (1976) 287. 

[3] N.D. Birrel and P.C.W. Davies, Quantum Fields in Curved Space, Cambridge Uni- 
versity Press (1982). 

[4] W. G. Unruh and R. M. Wald, Phys. Rev. D 29 (1984) 1047 

[5] P. G. Grove, Class. Quantum Grav. 3 (1986) 801 

[6] D. Raine, D. Sciama and P. Grove, Proc. R. Soc. A435 (1991) 

[7] W. G. Unruh, Phys. Rev. D 46 (1992) 3271 

[8] F. Hinterleitner, Ann. Phys. (N.Y.) 226 (1993) 165 

[9] S. Massar, R. Parentani and R. Brout, Class. Quantum Grav. 10 (1993) 385 

[10] J. Audretsch and R. Miiller, Phys. Rev. D 49 (1994) 4056; Phys. Rev D 49 (1994) 
6566; Phys. Rev A 50 (1994) 1755 



53 



[11] R. Brout, S. Massar, S. Popescu, R. Parentani and Ph. Spindel, " Quantum Source 
of the Back Reaction on a Classical Field' Preprint ULB-TH 93/16, UMH-MG 
93/03, (1993). 

[12] Y. Aharonov, D. Albert, A. Casher and L. Vaidman, Phys. Lett. A 124 199 (1987), 
Y. Aharonov and L. Vaidman, Phys. Rev. A 41 11 (1990). 

[13] J. von Newmann, Mathematical Foundations of Quantum Mechanics, Princeton 
University Press, Princeton (1955) 

[14] J.M. Bardeen, Phys. Rev. Letters 46 (1981) 382. 

[15] R. Parentani and T. Piran, Phys. Rev. Lett. 73 (1994) 2805 

[16] S. Massar, The semi classical back reaction to black hole evaporation preprint ULB- 
TH 94/19, |sr-qc/941 1039| 

[17] G. 't Hooft, Nucl. Phys. B 256 (1985) 727 

[18] T. Jacobson, Phys. Rev. D44 (1991) 1731, D48 (1993) 728. 

[19] F. Englert, S. Massar and R. Parentani, Class. Quantum. Grav. 11 (1994) 2919 

[20] L. Susskind, Phys. Rev. D 49 (1994) 6606 

[21] L. Susskind, "Some speculations about Black Hole Entropy in String Theory" 
Preprint RU-93-44, |hep-th/9309ll5| (1993). 

[22] L. Susskind, L. Thorlacius and J. Uglum, Phys. Rev. D 48 (1993) 3743 

[23] C. R. Stephens, G. 't Hooft and B. F. Whiting, Class. Quant. Grav. 11 (1994) 621 

[24] G. 't Hooft, Horizon Operator Approach to Black Hole Quantization preprint THU- 
94/02 (1994) |gr-qc/9402037| 

[25] F. Englert, "Operator weak values and black hole complementarity" , preprint ULB- 
TH 03/95, to be published in the annals of the Oskar Klein Centenary Symposium 

[26] K. Schoutens, H. Verlinde, E. Verlinde, "Black Hole Evaporation and Quantum 
Gravity" Preprint CERN-TH. 7142/94, PUPT-1441, (1994), |hep-th/ 94010811 . 

[27] P. Kraus and F. Wilczek, Nucl. Phys. B433 (1995) 403 

[28] D. G. Boulware, Annals of Physics 124 (1980) 169 

[29] R. Parentani, Class. Quantum Grav. 10 (1993) 1409. 

[30] R. Parentani and R. Brout, Int. J. Mod. Phys. Dl, 169 (1992). 



54 



[31] R. Brout, S. Massar, R. Parentani and Ph. Spindel, A Primer for Black Hole Quan- 
tum Physics (1995) ULB-TH 95/02, UMH-MG 95/01 and LPTENS 95/03 submitted 
to Phys. Rep. 

[32] R. Parentani, The Recoils of the Accelerated Atom and the Decoherence of its Fluxes, 
preprint (1995) LPTHENS 95/02 

[33] R. Wald, Commun. Math. Phys. 45 (1975) 9. 

[34] P. G. Grove in The Origin of Structure in the Universe edited by E. Gunzig and P. 
Nardone, Kluwer Academic Publishers (Netherlands) 1993 

[35] R. Wald, Commun. Math. Phys., 54,1 (1977), Phys. Rev. D, 17, 1477 (1978), 

[36] S. Massar, R. Parentani and R. Brout, Class. Quantum Grav. 10 (1993) 2431 

[37] S. Massar, Int. J. Mod. Phys. D 3 (1994) 237 



55 



Figure Captions 
Figure 1. 

The Minkowski coordinates t, z and U, V. The left (L) and right (R) Rindler quadrants. 
The Rindler coordinates r, p in R and the trajectory of a uniformly accelerated atom. 
Figure 2. 

The absolute value of the switch function /(r) given in eq. fl6"4|) for m = 2a and T = 3a _1 . 
r is given in units of a -1 . 
Figure 3. 

The mean Rindler energy density ((T vv (v)) therm.) emitted to order g 2 at thermal equi- 
librium is represented for m = 2a and T = 3a" 1 . v is given in units of a -1 and T vv in 
arbitrary units since the flux is proportional to the coupling g. One sees the vanishing 
of the flux in the steady regime and the positivity of the transients. In the Minkowski 
description they are enenhanced by the jacobian dV/dv to make the total Minkowski 
energy emitted positive. 
Figure 4. 

The conditional value {T vv ) e if the two level atom is initially in its ground state and ends 
up in its excited state. The parameters are the same as in Figure 2 and 3: m = 2a and 
T = 3a _1 . The v axis is given in units of a" 1 and T vv in units of a 2 . For U < U a , V > 0, 
(T vv ) e is complex and oscillates. The real part has a central positive bump which encodes 
that their is a rindleron carrying positive energy which will induce the transition of the 
atom. For V < 0, (T vv ) e is real and positive. It describes the partner of the rindleron 
which will be absorbed by the atom. The oscillations of (T VV (U < U a , V > 0)) e are such 
that the total Minkowski energy of the vacuum fluctuation vanishes. For U > U a , V > 0, 
(T vv ) e is positive and of order N m . In order to represent it we have had to change the 
vertical scale. 
Figure 5. 

A schematic picture of the energy fluxes (T m ) e . We have represented in dark grey the 
regions where (T vv ) e is 0(N m + 1) and in light grey the regions where it is 0(N m ). 
Figure 6. 

The local description of a vacuum fluctuation giving rise to a Hawking photon emitted 
around u = uq is represented in a Penrose diagram. The shaded areas correspond to 
the regions where T^(-u) is non vanishing, v — v$ is the trajectory of the collapsing 
spherically symmetric shell of massless matter. 
Figure 7. 

The same as in figure 6 drawn in Eddington-Finkelstein coordinates. 
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